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1.  Introduction 


We  present  a  new  approach  for  the  construction  of  wavelets  and  prewavelets  on  IR*  from 
multiresolution.  Our  method,  which  is  based  on  our  earlier  work  [BDR]  and  [BDR1],  uses  only 
properties  of  shift-invariant  spaces  and  orthogonal  projectors  from  £^(01*)  onto  these  spaces,  and 
requires  neither  decay  nor  stability  of  the  scaling  function.  Furthermore,  this  approach  allows  us 
to  derive  in  a  simple  way  previous  constructions  of  wavelets,  as  well  as  new  constructions,  and  to 
settle  completely  certain  basic  questions  about  mult'resolution. 

A  univariate  function  ip  £  £2(01) 1S  called  an  orthogonal  wavelet  if  its  normalized,  translated 
dilates  ipjtk  :=  2k/2ip(2k  •  —  j),  j,  k  £  5Z,  form  an  orthonormal  basis  for  .^(IR).  In  other  words,  this 
system  is  complete  and  satisfies  the  orthogonality  conditions 

(1.1)  f  ipj'kipj.'k,  =  6(j  -  j')6(k  -  fc'),  j,k,j',k'  £  7Z, 

with  (5  the  delta  function  on  2Z.  The  concept  of  prewavelet  is  somewhat  more  general  in  that 
it  requires  (1.1)  to  hold  only  when  k  /  k'  and  hence  the  functions  there  are  not  assumed  to  be 
ortiiogonal  at  a  fixed  dyadic  level  k.  In  particular,  ip{-  —  j ),  j  £  7h,  are  not  necessarily  orthogonal, 
and,  instead,  one  assumes  that  (ip(-  —  j))je%  forms  a  stable  basis  for  £,2(IR)  (see  the  end  of  this 
section  and  §2  for  the  definition  of  stability). 

On  HT*,  wavelet  and  prewavelet  bases  are  generated  by  the  translation  and  dilation  of  the 
elements  of  a  set  $  of  2d  -  1  functions  from  £,2(0t*).  We  shall  say  that  they  are  an  orthogonal 
wavelet  set  if  {rpj,k  ■=  2kd/2ip(2k  •  -  j)  :  ip  £  '£,j  £  7Ld,k  £  ZZ)  is  an  orthonormal  basis  for 
£,2(111*).  Analogously,  the  set  $  is  a  prewavelet  set  if  {i pjtk  :  ip  €  €  7Zd,k  €  2Z}  is  a  stable 

basis  for  F2(DId),  and  in  addition  we  have  orthogonality  between  levels: 


(1.2)  /  iphk<p},'k.  =0,  k  ±  k\  jj'£7Ld ,  IMG*- 

The  construction  of  orthogonal  wavelets  has  a  rich  history  described  in  the  monograph  of 
Meyer  [Me]  and  the  article  of  Daubechies  [Dl].  Prewavelets  have  been  stressed  only  in  recent 
years  beginning  with  Battle  [B] .  The  paper  of  Jia  and  Micchelli  [JM]  discusses  their  brief  history. 
Most  methods  used  for  the  construction  of  wavelets  are  based  on  the  notion  of  multiresolution  as 
introduced  by  Mallat  [Ma]  and  Meyer  (see  [Me]).  Muitiresolution,  which  we  now  describe,  will  also 
form  the  starting  point  for  our  constructions. 

We  say  that  a  space  S  of  functions  defined  on  IRrf  is  shift- invariant  if,  for  each  s  £  S,  the 
shifts,  s(-  -  j),  j  £  7Ld ,  of  $  are  also  in  S.  More  generally,  we  say  that  S  is  /i-shift-invariant  if  it  is 
closed  under  /iK1*- translations.  All  shift-invariant  spaces  considered  in  this  paper  are  assumed  to 
be  closed  subspaces  of  £,2(01*).  Important  examples  of  shift-invariant  spaces  are  those  generated 
by  a  finite  set  $  :=  {<pi of  functions  from  £,2(01*).  For  such  $,  we  define 

S*($) 

to  be  the  £,2(IR*)-closure  of  the  finite  linear  combinations  of  the  2~fc-shifts  of  the  functions  from 
$.  We  write 


S($)  :=  S°($). 


Av;„i Aafcixity  Code# 

|Ava41  and/or 
at  .  Special 


□  □ 


In  case  $  consists  of  a  single  element  <j> ,  we  write  Sk(</>)  (instead  of  <Sfc({<£}))  and  we  say  that  S(<p) 
is  a  principal  shift-invariant  space. 

Now  suppose  that  we  hold  in  hand  a  sequence  of  spaces  {«S*}fc€z,  with  Sk  a  2-fc-shift-invariant 
space  for  each  k  €  ZZ.  We  shall  say  that  {5*}  forms  a  multiresolution  if  the  following  conditions 
are  satisfied: 

(i)  C  Sk+1,  k  €  7L- 

(1.3)  (ii)  US k  =  L2{m.dy, 

(iii)  r\Sk  =  {o}. 

In  the  usual  definition  of  multiresolution  analysis  as  proposed  by  Mallat  [Ma],  it  is  also  assumed 
that  (a)  Sk  is  the  2fc-dilate  of  some  fixed  principal  shift-invariant  space  and  that  (b)  the  shifts 
of  4>  form  an  orthonormal  family.  We  do  not  assume  these  conditions  in  our  definition  (1.3)  in 
order  that  we  can  discuss  more  general  situations  that  are  covered  by  the  techniques  of  this  paper. 
However,  for  the  remainder  of  this  introduction,  in  order  to  keep  the  discussion  simple,  we  shall 
assume  (a),  i.e.,  that  Sk  is  of  the  form 

=  (s(2fc.) :  s  €  S(<f>)}, 

for  some  principal  shift-invariant  space  «S(0).  Equivalently,  each  Sk  is  generated  by  the  2_fc-shifts 
of  the  dilated  function  4>(2k-).  In  this  case,  the  condition  (1.3)(i)  is  already  implied  by 

(1.4)  S°CSl. 

Previous  constructions  of  wavelets  assume  that  <p  has  /^(IR^-stable  shifts,  a  notion  which  we 
now  introduce.  For  a  collection  F  C  Z-2(IRd),  we  say  that  F  is  a  stable  basis  (for  the  space  it 
generates)  if  there  exist  positive  constants  C\{F)  and  C2(F)  such  that,  for  any  finitely  supported 
a  :=  {af) j£F, 

(1-5)  C'lC'OlHktF)  <  II  Z  “//II  <  ct(f )IHI«,(f). 

/€F 

In  the  context  of  wavelets,  the  family  F  is  taken  to  be  the  2-fcZd-shifts  of  some  finite  set  4>,  with 
the  totality  of  shifts  taken  over  all  k  6  TL  or  for  some  fixed  k.  Discussions  of  the  stability  question, 
including  earlier  references,  can  be  found  in  [JM]  (also  for  norms  other  than  the  2-norm)  and  in 
[BDR1].  Because  the  finitely  supported  sequences  are  dense  in  l2(7Ld),  the  Instability  of  the  shifts 
of  a  function  <f>  implies  that  the  map 

t2{7Ld)  -  S{4>) :  a  ^  £  4>(- -  j)a(j) 

jZTL* 

is  well-defined  and  induces  a  Hilbert  space  isomorphism  between  l2(7Ld)  and  S{4>). 

Thus,  under  the  stability  assumption,  (1.4)  is  equivalent  to  having  a  refinement  equation 

(1.6)  <f>(x)  =  ^(2l  -  iH?) 
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hold  for  some  sequence  a  £  C2(7Zd),  or,  equivalently,  to  have 

$  =  M(-/2) 

for  some  A  £  L2(2TTd),  called  the  (refinement)  mask. 

Given  a  shift-invariant  space  S  :=  <S°  whose  2fc-dilates  Sk ,  k  £  7L,  satisfy  (1.3)(i),  we  define 
the  wavelet  space  W  as  the  orthogonal  complement  of  5°  in  S 1 : 

W  :=  S1  QS°. 


It  follows  that  Wk  :=  Sk+l  ©  Sk ,  k  £  TL,  is  the  2fc-dilate  of  W.  The  spaces  Wk ,  k  £  7L,  are 
mutually  orthogonal. 

One  can  equivalently  define  W  by  projections.  If  P  Ps  is  the  orthogonal  projector  from 
/^(HU1)  onto  S,  then  W  =  {s  —  Ps  :  s  £  «Sl } .  If  (1.3)(ii)  and  (iii)  are  also  satisfied,  then  one  has 
the  orthogonal  decomposition 

(1-7)  I2(IRd)=®IFfc 

fc€2Z 

since,  for  each  /  6  LifJR?),  we  have 


/=  E(pfc/-pfc-i/).  with  Pkf-Pk-xf  ewk-\k£7L, 

fce  2Z 


and  Pk  the  Z2(IR<<)  projector  onto  Sk  (which  is  obtained  from  P  by  dilation).  Indeed,  the  condition 
(1.3)  (ii)  implies  that  limfc_00  Pkf  =  /,  and  (1.3)(iii)  implies  that  ]im k-. -oo  Pkf  =  0.  Note  that 
Pk  -  Pk- 1  is  the  orthogonal  projector  of  onto  W*-1. 

One  obtains  wavelets  and  prewavelets  from  multiresolution  analysis  by  finding  generators  for 
the  space  W.  For  example,  in  the  univariate  case,  Mallat  [Ma]  begins  with  a  function  <j)  £  I2(IR) 
which  has  orthonormal  shifts  and  satisfies  (1.3)  (with  Sk  :=  Sk((j>))  and  shows  that  W  is  a  principal 
shift-invariant  space  W  =  with  xp  an  orthogonal  wavelet.  One  can  also  apply  the  Mallat 

construction  to  a  function  4>  whose  shifts  are  only  T2(Dl)-stable  by  first  orthonormalizing  these 
shifts. 


Unfortunately,  if  <f>  is  of  compact  support  (and  its  shifts  are  not  orthonormal)  then  the  orthog¬ 
onal  wavelet  ip  will  generally  not  have  compact  support.  This  motivated  the  study  of  prewavelets. 
We  obtain  prewavelets  xp  by  finding  generators  of  W  whose  shifts  form  an  T2(IR)-stable  basis  for  W 
(but  not  necessarily  an  orthonormal  system).  Chui  and  Wang  [CW]  and  Micchelli  [Mi]  have  shown 
in  the  univariate  case  that  if  <p  has  compact  support  and  I2(1R)-Stable  shifts  and  (1.3)  is  satisfied 
(again  with  Sk  :=  Sk(<p ))  then  there  is  a  compactly  supported  prewavelet  x}x  which  generates  W . 
Chui  and  Wang  even  characterize  the  xp  £  W  of  minimal  support  (in  a  sense  to  be  made  clear  in 
§5)  which  generates  W.  We  shall  give  a  simple  derivation  of  (a  slightly  stronger  version  of)  these 
facts  in  §5. 

In  the  multivariate  case,  the  construction  of  orthogonal  wavelet  and  prewavelet  sets  is  by  far 
more  involved.  Micchelli  [Mi]  and  Jia  and  Micchelli  [JM]  have  studied  multiresolution  in  the  case 
the  function  <p  has  L2(IRrf)-stable  shifts  and  satisfies  two  regularity  conditions.  The  first  of  these  is 
that  the  periodization 


(1.8) 


£  w  -j)i 

j€Z‘i 
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of  \4>\  is  in  L2{T£d).  (Note  that  this  requirement  is  satisfied  if  <j>  has  suitable  decay  at  oc).  Secondly, 
they  require  that  (f)  satisfies  the  refinement  equation  (1.6)  with  the  coefficient  sequence  a  in  t\{TLd). 

In  contrast  with  the  present  literature,  we  only  need  to  assume  here  that  the  function  <f)  is  in 
£2(IRd),  satisfies  the  refinement  condition  (1.3)(i),  (with  <S°  :=  S(<f>),  and  Sk  the  2fc-dilate  of  S° ) 
and  that  its  Fourier  transform  4>  satisfies 


(1.9)  supp$=IRd. 

Here  and  later,  the  support  of  an  ^(K^Munction  /  is  defined  only  modulo  a  null-set  as  {x  :  }(x)  / 
0}  and  the  Fourier  transform  of  a  function  /  €  Ti(IRd)  is  defined  by 


(1.10) 


e-y  ft 


where,  here  and  throughout, 

c$  :  x  i — *  e'x'6 

is  the  complex  exponential  with  frequency  0  6  IIId.  (We  shall  use  without  further  mention  basic  facts 
from  Fourier  analysis  including  the  fact  that  the  Fourier  transform  has  an  extension  to  Z2(flld).) 
In  particular,  our  analysis  applies  whenever  4>  has  compact  support  since  then  <f>  is  'nalytic  and 
its  zero  set  is  of  measure  zero  (unless  (f>  =  0).  We  note  that  we  shall  not  need  to  assume  that  <j> 
has  l2(Dfrf)-stable  shifts,  nor  impose  any  decay  conditions,  nor  any  conditions  on  the  refinement 
coefficients  a.  In  fact,  we  do  not  even  need  to  assume  the  refinement  condition  in  the  form  (1.6), 
only  in  the  original  form  (1.3)(i). 

Under  the  above  assumptions,  we  show  in  this  paper  that  the  conditions  (1.3)(ii)  and  (1.3)(iii) 
of  multiresolution  automatically  hold.  Further,  our  derivation  of  (1.3)(iii)  from  (1.3)(i)  does  not 
make  use  of  (1.9).  We  even  provide  a  characterization  of  property  (1.3)(ii)  for  the  case  when  (1.9) 
fails  to  hold.  We  also  show  that  (1.3)(ii)  and  (iii)  automatically  hold  whenever  <p  is  of  compact 
support.  Details  can  be  found  in  §4.  Previous  results  on  the  matter  (cf.  e.g.,  [JM]  and  [So])  were 
derived  under  the  stability  assumption  and  under  suitable  decay  conditions. 

The  main  goal  of  multiresolution  is  to  construct  a  set  of  $  of  2d  -  1  functions  which  generate 
the  wavelet  space  XV  (i.e.,  W  =  «S(^))  and  have  other  prescribed  properties.  We  shall  index  the 
elements  in  by  the  set  V1  :=  V  \  {0},  with  V  the  set  of  vertices  of  the  cube  [0  . .  l/2]rf.  A  major 
advantage  of  our  approach  is  that  it  is  almost  trivial  to  find  generating  sets  'J'  for  W.  Once  one 
such  set  ^  is  found,  we  can  then  find  (all)  other  generating  sets  by  simple  operations  on  the  Fourier 
transforms  of  the  elements  of  9. 

Two  particularly  interesting  generating  sets  which  are  obtained  by  our  construction  are  dis¬ 
cussed  in  §3.  First  of  all,  we  show  that  (1.9)  implies  that  W  always  possesses  a  generating  set  $ 
which  provides  an  orthonormal  basis  for  W,  i.e.,  an  orthogonal  wavelet  set.  Secondly,  under  slightly 
more  restrictive  assumptions  on  <f>,  we  show  that  there  is  a  function  w  6  L 2(IRd)  whose  half-shifts 
tu(-  +  v),  v  €  V1 ,  form  a  generating  set  for  W.  Special  cases  of  this  latter  result  have  been  proved 
in  [MRU]  and  [Mil],  see  also  [LM]. 
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A  more  delicate  problem  is  the  construction  of  multivariate  wavelets  and  prewavelets  which 
have  compact  support.  If  the  function  4>  of  multiresolution  has  compact  support,  it  is  quite  easy  to 
find  generating  sets  for  W  whose  elements  are  compactly  supported.  On  the  other  hand,  if  <f>  has 
Z2(IR‘i)-stable  shifts,  we  would  like  the  shifts  of  the  functions  in  S'  to  form  an  Z2(IRd)-stable  basis 
for  W .  While  it  has  been  shown  by  Meyer  [Me;  Chapter  III,  §6],  and  Jia  and  Micchelli  [JM1]  that 
such  generating  sets  always  exist,  their  proofs  are  not  constructive.  On  the  other  hand,  several 
authors,  including  Riemenschneider  and  Shen  [RSI],  Chui,  Stockier,  and  Ward  [CSW],  Lorentz 
and  Madych  [LM],  and  Stockier  [So],  have  given  constructions  of  prewavelet  sets  S'  under  various 
conditions  on  cj)  and  in  some  cases  with  restrictions  on  the  space  dimension  d.  We  shall  discuss 
this  question  in  §7  where  we  shall  use  our  characterizations  of  the  wavelet  space  W  to  recover  and 
slightly  improve  some  of  these  constructions. 

A  particularly  interesting  application  of  wavelet  constructions  is  to  functions  (f>  which  are 
B-splines  or  box  splines.  In  this  regard,  we  obtain  the  compactly  supported  univariate  spline 
prewavelets  of  Chui  and  Wang  [CWl]  and  derive  various  orthogonal  wavelets  and  prewavelets 
obtained  from  box  splines. 

As  we  have  already  noted,  our  construction  of  wavelets  and  prewavelets  is  based  on  our  earlier 
results  on  the  structure  of  shift-invariant  spaces.  We  use  two  facts  repeatedly.  The  first  is  an  exact 
description  of  finitely  generated  shift-invariant  spaces.  For  example,  we  have  shown  in  [BDR]  that 
the  principal  shift-invariant  space  <S(d>)  is  described  by  its  Fourier  transforms: 

(1.11)  <S(d>)  =  {T(j>  €  Z2(lR<i)  :  T  is  27r-periodic). 

Here  and  later,  for  a  set  of  functions  F,  we  define  F  :=  {/  :  /  €  F}  to  be  the  set  of  its  Fourier 
transforms.  A  similar  characterization  holds  for  a  finitely  generated  space  (see  §2). 

In  the  case  that  d>  has  Z2(Dld)-stable  shifts,  (1.11)  is  well-known,  and  the  functions  r  must  be 
in  L2{TTd),  with 

TT* 

the  d-dimensional  torus.  Assuming  suppd>  =  IR  ,  we  have  shown  in  [BDR1]  that  there  is  always 
a  function  d>„  which  generates  S(4>),  (i.e.  5(d>.)  =  <S(d>)),  whose  shifts  are  Z2(lRd)-stable;  in  fact 
they  can  be  taken  to  be  orthonormal. 

The  second  result  which  we  frequently  employ  is  the  explicit  formula  (2.11)  of  the  next  section 
for  the  orthogonal  projector  P  =  P#  from  Z2(lR‘i)  onto  the  principal  shift-invariant  space  S(4>). 

With  these  results  in  mind,  our  construction  proceeds  as  follows.  We  show  that  if  d>  is  in 
Z2(IRt/)  with  suppd>  =  HT*  and  if  the  space  sequence  (tS*(<£))fc62Z  sat>sfies  ( 1.3)(i),  then  we  can  give 
an  alternate  description  for  Sl: 

(1.12)  5!=5($) 

where  $  :=  (<£(•  +  v))v€V-  It  follows  that 

wv  :=  <f>(-  +  v)  -  F*(<£(-  +  v)),  reF', 
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are  a  set  of  generators  for  W .  We  then  use  our  characterization  of  finitely  generated  shift-invariant 
spaces  to  obtain  other  generators  with  more  favorable  properties.  Because  of  our  description  of 
the  projector  P#,  all  these  generators  are  described  in  a  concrete  way  in  terms  of  their  Fourier 
transforms.  In  this  way,  we  are  able  to  construct  an  orthonormal  basis  for  the  wavelet  space  without 
using  either  the  refinement  equation  or  the  mask,  and  further,  without  making  any  assumption  on 
the  stability  or  orthogonality  of  the  shifts  of  <f>.  Under  further  assumptions  (e.g.,  when  (f>  is  compactly 
supported),  generators  of  the  wavelet  space  that  can  be  written  as  a  finite  linear  combinations  of 
<t>(2  •  -j),  j  €  7Ld,  are  obtained. 

As  already  alluded  to  in  our  definition  (1.3),  we  shall  actually  work  in  the  more  general  setting 
of  non-stationary  wavelets  in  this  paper,  which  means  that  our  spaces  Sk ,  while  still  being 
assumed  to  be  generated  by  the  2-k-shifts  of  some  function  </>*,  will  not  be  assumed  to  be  the 
dilate  of  S°  or  of  any  other  space  in  the  sequence  {<SJ}.  It  turns  out  that  this  generalization  can 
be  handled  at  no  additional  cost  and  leads  to  interesting  bases  for  For  example,  in  §6, 

we  discuss  such  an  example  based  on  exponential  splines,  and  in  §8  w’e  discuss  their  multivariate 
analog,  the  exponential  box  splines. 

An  outline  of  the  present  paper  is  as  follows.  In  §2,  we  review  and  extend  results  from  our 
earlier  work  which  will  be  needed  in  the  sequel.  In  §3,  we  describe  generating  sets  and  bases  for 
wavelet  spaces  W .  In  §-!,  we  analyze  conditions  (1.3)(ii)  and  (iii)  of  multiresolution.  In  §5,  we  apply 
our  constructions  to  derive  univariate  wavelets  and  prewavelets  with  various  desirable  properties. 
In  §G.  we  discuss  exponential  B-splines  as  wavelets.  In  §7,  we  consider  the  construction  of  wavelets 
in  the  multivariate  case.  We  conclude  with  a  brief  discussion  in  §8  of  exponential  box  splines  as 
wavelets,  and  describe  stable  bases  for  their  associated  wavelet  spaces. 

2.  Shift-invariant  spaces 

Our  analysis  will  be  based  on  the  structure  of  shift-invariant  spaces  given  in  our  earlier  work 
[BDR]  and  [BDR1].  In  this  section,  we  review  some  of  these  facts  and  develop  them  somewhat 
further  in  directions  pertinent  to  the  construction  of  wavelets.  We  have  already  mentioned  in  (1.11) 
a  characterization  of  S  =  S((f>)  in  terms  of  Fourier  transforms.  A  similar  characterization  of  the 
space  *Sfc(d>),  generated  by  the  2~fc-shifts  of  (f> ,  easily  follows  from  (1.11)  by  dilation: 

(2.1)  Sk(<f>)  =  {t4>  €  £2(IRd)  :  t  is  2Ar+17r-periodic}. 

In  the  context  of  the  principal  shift-invariant  space  S(<f>),  it  will  be  important  to  know  whether 
some  given  function  /  €  S(4>)  generates  this  space,  i.e.,  whether  S(f)  =  With  the  aid  of 

(1.11),  we  obtain  the  following  simple  answer  [BDRl]  to  this  problem: 

Corollary  2.2.  Let  S(4>)  be  a  principal  shift-invariant  space,  and  let  J  £  S(d>).  Then  f  generates 
S(4>)  if  and  only  if  supp  /  D  supp0. 

Proof.  If  4>  €  S(f),  there  exists,  by  (1.11),  a  27r-periodic  r  such  that  <f>  =  r/,  and  hence 

.•** 

supp  /  D  supp  <{>. 
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For  the  converse,  we  assume  that  supp  /  3  supp  4>  and  want  to  show  that  S(f)  =  S(<p).  Since  we 
assume  that  /  £  S(<p),  then,  by  (1.11),  there  exists  r  such  that  /  =  a.e.  Defining  r'  to  be  1/r  on 
suppr  and  0  elsewhere,  we  obtain  that  a.e.  4>  =  r  /  on  suppr,  but  since  suppr  3  supp  /  j  suppd> 
(the  last  inclusion  by  assumption),  the  equality  holds  everywhere.  By  (1.11),  we  conclude  that 
<j>  £  S(f),  and  hence  S(f)  =  <$(<£).  4 

The  above  description  of  principal  shift-invariant  spaces  in  terms  of  their  Fourier  transforms 
can  be  generalized  to  the  finitely  generated  shift-invariant  space  5($)  as  follows  (cf.  Theorem  1.6 
of  [BDR1]): 

(2.3)  5(4?)  =  6  I2(II ld)  :  r $  is  2x-periodic,  <j>  £  $}. 

From  the  description  (1.11)  of  the  principal  shift-invariant  space  S(<f>),  we  see  that  the  Fourier 
transform  of  a  function  s  £  S(4>)  is  determined  by  its  values  on  TTd  (at  least  when  suppd>  =  IR^). 
It  is  possible  to  factor  out  this  redundancy  with  the  aid  of  the  bracket  product  which  is  defined 
for  /, g  €  Z.2(IR'i)  by  the  formula: 

(2.4)  [/, 3]  :=(/?)"=  Y.  /(•+«•  +  «• 

pZ'lxTL'1 

Note  that  [/, <7]  is  in  Li(TTd)  whenever  f,g  £  /^(IR1*).  Also,  by  the  Cauchy-Schwarz  inequality, 

(2.5)  \[f,g]\'  <  [/>/]($,  $]> 

with  the  right  side  finite  a.e.  The  importance  of  (2.4)  lies  in  part  in  the  following  formula,  valid 
for  f,g£  L2(TRd): 

(2-6)  Lj(x-mz)ix^^L€->n=wfLe-Ar'% 

which  shows  that  the  inner  product  of  /(•  —  j)  with  g  is  the  j-th  Fourier  coefficient  of  [/,</]. 

It  is  easy  to  derive  the  following  three  elementary  properties  of  the  bracket  product.  The  first 
two  follow  from  (2.6),  while  the  third  one  follows  directly  from  the  definition  (2.4). 

Lemma  2.7.  If  f,g  £  Z2( IR^),  then  the  shifts  of  f  are  orthogonal  to  the  shifts  of  g  if  and  only  if 

[/.?]  =  0. 

Lemma  2.8.  If  f,g  £  Z,2(fll'1)  are  compactly  supported,  then  [f,g]  is  a  trigonometric  polynomial. 

Lemma  2.9.  If  f,g  £  L2(IRd)  and  r  has  period  2n,  then  {rf,g}  =  r[f,g j  =  [/,  fg). 

The  bracket  product  also  appears  naturally  in  the  computation  of  the  norms  of  elements 
s  £  S(4>).  By  (1.11),  s  =  r<{>  and 

(2-10)  (2*')‘1/,2||s||L3(iv')  =  IIS11  J)  = 
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There  is  a  simple  description  for  the  orthogonal  projector  P  :=  from  L2(IRd)  onto  the 
principal  shift-invariant  space  S(4>).  For  each  /  €  Z.2(IRd),  P<pf  is  the  best  /^(IR^j-approximation 
to  /  from  S(<t> )  and  is  characterized  by  the  orthogonality  of  the  error  f  -  P<pf  to  S(<p).  It  was 
shown  in  [BDR]  that  P#  is  described  by 

(2.11) 

[<p,  <p\ 

where  we  use  the  convention  (throughout  this  paper)  that  0  times  any  extended  number  is  0;  in 
particular  0/0  is  defined  to  be  0.  (We  note  that,  by  the  definition  of  the  bracket  product,  <p  vanishes 
whenever  [<£,  does.)  There  is  a  similar  formula  (which  we  shall  not  need)  in  the  case  S  is  finitely 
generated  (see  [BDRl]). 

There  are  several  interesting  points  to  be  made  about  bracket  products  and  the  projector  P#. 

^  ^  i 

First  of  all,  from  (2.6),  it  follows  that  <£  has  orthonormal  shifts  if  and  only  if  [</>,  </>]  =  1  a.e.  on  TT  . 
In  this  case,  formula  (2.11)  is  the  (Fourier  transform  of  the)  usual  one  for  projecting  onto  S(<j>).  In 
the  case  that  <p  does  not  have  orthonormal  shifts,  but  [<£,  d>]  ^  0  a.e.,  the  function  <}>,  with  Fourier 
transform 


(2.12) 


is  in  S(4>),  has  orthonormal  shifts,  and  generates  «S(d>),  i.e. ,  S(<f>.)  =  S(4>). 

Since  the  square  root  of  the  bracket  product  appears  very  frequently,  we  introduce  the  following 
notation: 


(2.13)  £:=[£,$]*/2  =  (  Y,  \&'  +  P)\ 2),/2’ 

P&nTZ* 


with  the  sum  taken,  offhand,  pointwise,  hence  defined  for  any  <f>  on  Dld  if  we  allow  it  to  take 
the  value  oo.  From  (2.6),  we  conclude  that  the  map  Z.2(IRd)  — ♦  Z.2(Trd)  :  <f>  >-*  cf>  is  non-linear, 
norm-preserving,  and  onto: 

Lemma  2.14.  The  function  <t>  is  in  Lt( Dld)  if  and  only  if  €  Z,2(TTd).  Moreover,  ||d>||z,j(R'*)  = 
Turning  back  to  the  orthogonal  projection,  we  can  write  it  in  the  form 


(2.15)  P*f  =  [/,  £]</>,  £:=4"’ 

<P 

and  check  that  =  1  on  supp<£,  and  therefore,  by  the  above  lemma,  Ji  €  Z-2(IRd)  if  and  only 
if  l/</  6  L2(supp<£).  In  such  a  case,  by  (2.5),  [/,/ij  e  Zi(TTd),  and  we  can  formally  write  the 
orthogonal  projection  P^f  in  the  form 

(2.16)  P*/=  Y  ^ 

jai* 
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with  [/,p]v(j)  the  iourier  coefficient  of  [/,£]. 

A  special  case  of  the  above  occurs  when  4>  has  ^(IR^-stable  shifts,  i.e.,  when  (<£(•  +  j))je7ld 
forms  a  stable  basis  for  As  will  be  explained  below,  this  stability  is  equivalent  to  having 

Cl  <  4>  <  C2  a.e.  on  TT^  for  constants  C\,Ci  >  0.  In  this  case,  the  formed  expansion  (2.16) 
cor  .erges  (unconditionally)  and  hence  the  orthogonal  projection  takes  the  explicit  form 


(2.17) 


p,i  =  £  <«• 

jell4 


+  j)nj(j), 


m/0)“  / 

Jnd 


f(x)v{x  +  j)dx. 


A  similar  analysis  applies  to  the  structure  of  the  spaces  £*(<£)  generated  by  the  2~fc-shifts  of 
the  function  <f>.  We  shall  only  need  the  case  £'(<£)  in  the  sequel.  In  this  case,  the  bracket  product 
is  replaced  by  the  double-bracket  product 

(2-18)  =  Y.  K  +  PM'  +  P)* 

064  w7ld 


which  is  a  4ff-periodic  function.  The  role  of  /  is  played  by 

(2.19)  /:=(  Y  l/(-  +  /3)|2)1/2,  fei2(md). 

P&inTLd 


Note  that  /  is  a  non-negative  Iff- periodic  function,  and 

(2.20)  J(xf  =  Y  /(J+/?)2- 

/3€4jtV' 

with  V  the  vertices  of  the  cube  [0  . .  1  /2]d. 

In  particular,  <t>  characterizes  stability  and  orthogonality  of  the  half-sKihs  of  4>.  Namely,  the 
orthogonality  of  the  half-shifts  is  equivalent  to  <f)  =  2~dC  a.e.,  and  the  stability  is  equivalent  to  the 
boundedness  a.e.  of  <*>  and  l/<£. 

We  next  describe  in  more  detail  the  structure  of  the  finitely  generated  space  S  =  <S(^);  a  more 
complete  discussion  can  be  found  in  [BDRlj.  First,  for  s  €  5,  the  representation  (2.3)  for  s  is  local 
in  the  sense  that  we  can  independently  assign  the  values  r^x)  for  x  6  TT^  and  </>€$.  The  choice 
of  r0(x)  determines  the  value  of  s  at  all  points  in  x  -f  2x77d .  This  means  that  the  structure  of  S  is 
determined  by  the  vectors  (^(x  -f  fl))t3£2x7Z.di  4>  €  $,  x  €  TT^.  For  example,  for  any  fixed  x,  these 
vectors  are  linearly  independent  if  and  only  if  the  associated  Gramian  matrix 

(2.21)  G<*>  ==(ft  «),,,» 


has  non-zero  determinant  at  x.  In  particular,  if  detG($)  =  0  on  a  set  of  positive  measure  in  'IT<i, 
then  the  representation  (2.3)  is  not  unique. 
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We  say  that  the  set  of  generators  $  provides  a  basis  for  S  (that  is,  their  shifts  are  a  basis)  if 
the  representation  (2.3)  is  unique  for  each  s  6  S,  or,  equivalently,  if  detG($)  is  nonzero  a.e.  All 
bases  for  S  have  the  same  number  of  elements.  We  note,  [BDR1],  that  not  every  finitely  generated 
shift-invariant  space  S  contains  a  basis.  We  also  note  that  G($)  is  a  non-negative  definite  matrix, 
hence  $  is  a  basis  if  and  only  if  det(?(<£)  >  0  a.e.  on  TTd. 

We  further  say  that  a  set  of  generators  $  provides  an  1,2-stable  basis  for  <S($)  if  each  s  £  S 
has  a  unique  representation 

(2.22)  3  =  S  <K‘  ~  j)ciM 

<t>e*  je7id 


and  the  coefficients  satisfy 

(2-23)  C\  lcj,<*>(s)|2 

<t>e*  jen4  *€*  jeTi* 

with  absolute  constants  Ci,C2  >  0.  Any  L2-stable  basis  is  obviously  a  basis.  One  easily  checks 
that  the  present  definition  coincides  with  the  (seemingly  weaker)  one  given  in  the  introduction. 

We  recall  from  [BDRl]  that  the  finite  generating  set  $  for  S  is  an  Z,2-stable  basis  for  5  if  and 
only  if,  for  some  matrix  norm  ||  •  ||  (and  hence  all  matrix  norms), 

(2.24)  IKWII,  \\G(*)-'\\  €  Ioo(TTd). 

In  particular,  this  is  the  case  only  if 

(2.25)  C\  <  det  G($)  <  C2,  a.e.  on  TT1*, 
for  absolute  constants  C\ ,  C2  >  0. 

In  the  case  that  the  entries  of  G($)  are  continuous,  the  ‘a.e.’  in  (2.25)  can  be  removed,  and 
more  importantly,  (2.25)  becomes  equivalent  to  the  Z,2-stability  of  $.  Furthermore,  in  this  case  the 
right  inequality  of  (2.25)  trivially  holds,  and  thus,  due  to  the  continuity  of  detG($),  stability  is 
equivalent  to  the  condition 

detC($)(x)  >  0,  Vi€Trd. 

This  latter  characterization  of  stability  was  obtained  by  Jia  and  Micchelli  [JM]  under  slightly 
stronger  assumptions  (which  imply  the  continuity  assumption). 

With  the  notions  of  basis  and  I2-stable  basis  in  hand,  the  following  theorem  shows  how  from 
one  $  which  provides  a  (stable)  basis  for  S  we  can  obtain  other  sets  with  the  same  property.  In  this 
theorem,  T  :=  (TV^Vetf^e*  denotes  a  square  matrix  whose  entries  are  2x-periodic  measurable 
functions. 

Theorem  2.26.  Let  $  provide  a  basis  for  S  =  <S($).  For  any  set  of  functions  from  S($),  we 
have: 

(i)  $  provides  a  basis  for  S  if  and  only  if  ’k  =  T$  for  some  T  which  is  nonsingular  a.e.; 

(ii)  'F  provides  a  basis  for  S  if  and  only  if  it  generates  S  and  =  #$; 
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(iii)  $  provides  a  basis  for  S  if  and  only  if  and  detG('f)  ^  0  a.e. 

(iv)  <5  provides  an  L^-stable  basis  for  S  if  $  does  and  $  =  T$  with  ||T||,  IjT-1 1|  in  L00(lTd). 

The  above  is  easily  proved  by  noticing  the  effect  of  the  transformation  T  on  the  Gramian: 

(2.27)  G((T$)V)  =  TG($)Tm, 

with  T*  the  conjugate  transpose  of  T;  see  [BDR1;  Cor. 3. 31]  for  more  details. 

A  special  case  of  the  above  occurs  when  T  is  a  diagonal  matrix.  In  this  case,  provides  a  new 
basis  if  and  only  if  ^  C  L 2(lRd)  and  the  2?r-periodic  functions  {r^}  are  different  from  zero  a.e. 
on  TT^.  Also,  $  is  an  I^-stable  basis  if  $  is,  and  the  and  their  reciprocals  are  in  Z/cc,('ITd). 

If  $  provides  a  basis  for  S  :=  <S($),  then  $  can  be  orthonormalized  by  a  Gram-Schmidt 
orthogonalization.  We  summarize  this  fact  in  the  following  theorem  whose  proof  is  left  to  the 
reader  (and  can  be  found  in  [BDR1]). 

Theorem  2.28.  Let  $  provide  a  basis  for  S  :=  *S($). 

(i)  Then,  there  is  a  set  4>*  of  generators  for  S  that  provides  an  orthonormal  basis  for  S. 

( ii )  If  the  functions  in  $  have  compact  support,  then  there  is  a  set  . .  .,<£*}  of  compactly 

supported  functions  which  give  the  orthogonal  decomposition: 


(2.29)  S  =  S(<i>\)®---®S(<f>mn). 

The  half-shift-invariant  space  Sx(<fi),  generated  by  <f>,  is  identical  with  £($),  $  :=  (</>(•  +  u)) v€y, 
and  V  the  vertices  of  the  cube  [0  . .  l/2]d,  as  before.  Clearly,  orthonormality  or  /^-stability  of  the 
half-shifts  of  <j>  is  the  same  as  orthonormality  or  /^-stability  of  the  full  shifts  of  4>.  Thus,  there 

must  be  a  relation  between  G($)  and  <f>  which  we  shall  now  derive. 

Given  a  47r-periodic  function  /,  the  functions 

(2.30)  <?„(/):=  Y  ev(.  +  p)f(-  +  p),  v  €  2Zd/2 

are  2x-periodic.  If  A  has  been  obtained  by  choosing  exactly  one  point  from  each  of  the  cosets 
t;  -f  TLd ,  v  €  V ,  then 

(2.31)  /  =  X)  *—<?«(/)/ 2d 

v€  A 


is  a  decomposition  of  /  into  its  2^-periodic  components  Qv{f)/ 2d,  v  £  A. 

Since  7Ld / 2  is  the  disjoint  sum  V  +  and  ev  is  47r-periodic  for  v  €  7Ld / 2,  we  find  that,  for 
v,u$V, 


(2.32) 


[e„<£,  eu<£]  =  Y  ev(.  +  /?)e_u(- +  /?)!<£(•  +  /?)  |2 

=  Y  ev-u(- +  +  v)2  =  Qv-u(d>2)- 


With  this,  we  can  easily  compute  the  eigenvalues  and  eigenvectors  of  G($). 
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Lemma  2.33.  For  each  n  €  4xV  and  x  €  TT^,  the  number  2d4>(x+n)2  is  an  eigenvalue  ofG($)(x), 
with  eigenvector  au  :=  ( ev(x  +  n))veV- 

Proof.  In  view  of  (2.32),  the  u-th  component  of  G($)aM  is 

53  Qv-uG>2)eu( ■  +  fi)  =  ev( ■  +  /i)  53  +  (*)Qv-uG>2)  =  2de„(-  +  n)4>(-  +  n)2 

ugv  uev 

where,  in  the  last  equality,  we  used  (2.31)  as  well  as  the  27r-periodicity  of  Qv-u{4>2)-  ♦ 

Corollary  2.34.  If  supp<£  =  IRd,  then  the  set  $  :=  (<£(•  +  u))u€V  is  a  basis  for  S($). 

Proof.  From  Lemma  2.33,  detG($)  =  const  $('  +  A1)2-  From  our  assumption  on 

the  support  of  4>,  it  follows  that  <£,  hence  detG($),  does  not  vanish  a.e.,  hence  $  is  a  basis  for 

S($).  ♦ 

3.  Generators  for  the  wavelet  space 

We  give  in  this  section  various  descriptions  for  wavelet  spaces  and  their  generators  and  bases. 
As  mentioned  earlier,  we  develop  this  analysis  in  the  following  more  general  framework  than  usually 
considered  in  the  multiresolution  construction  of  wavelets:  We  suppose  that  <f>  and  tj  are  functions 
in  L2(JR.d)  with  the  property 

(3.1)  S(4>)CS1(T}). 

As  before,  S(4>)  denotes  the  principal  shift-invariant  space  generated  by  <j>,  and  5J( tj)  denotes  the 
space  generated  by  the  half-shifts  of  tj.  Thus,  Sl(r])  is  also  the  2-dilate  of  S(t](-/2)). 

We  define  the  wavelet  space  W  as  the  orthogonal  complement  of  S(<f>)  in  S1(tj): 

W:=  s'wesit). 

The  purpose  of  this  section  is  to  find  generating  sets  and  bases  for  W.  The  analysis  of  this  section 
also  applies  to  the  other  wavelet  spaces 

(3.2)  Wk  :=  Sfc+1  0  Sk 

(which  might  be  generated  by  some  other  L2 -functions)  after  a  suitable  dilation.  In  the  case  of 
stationary  decompositions  usually  considered  in  wavelet  constructions,  we  have  rj  =  4>{ 2-),  and 
5'(i/)  is  the  2-dilate  of  «S(d>).  However,  (3.1)  is  much  more  general  since  we  can  begin  with  any 
Tj  e  £,2(1^)  and  take  for  <p  any  element  in  Sx{ij).  Indeed,  S{4>)  C  Sx(rj),  because  (77)  is  invariant 
under  shifting  by  half-integers,  a  fortiori  by  integers. 

Note  that  we  are  not  assuming  that  4>  is  rcfmable  in  terms  of  tj,  i.e.,  we  are  not  assuming  that 
<(>  can  be  written  in  the  form 

(3.3)  <t>=  53  *?(• -£>(■?') 
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for  some  sequence  a  (with  convergence  in  the  sense  of  Z^IR-'*))-  However,  since  S{4>)  C  (S1^),  we 
do  have 

(3.4)  4>=  Arj 

for  some  47r-periodic  function  A.  We  call  A  the  refinement  mask. 

We  do  assume  that 

(3.5)  supp  rj  =  supp  <f>. 

Although  this  assumption  is  not  essential  for  our  wavelet  constructions,  it  significantly  simplifies 
the  underlying  analysis.  In  case  <f>  and  77  are  compactly  supported,  as  is  the  case  in  almost  all 
wavelet  constructions,  (3.5)  automatically  holds,  since  then  supp  <f)  =  supp  rj  —  IR*. 

We  shall  now  describe  our  first  set  of  generators  for  W;  other  generating  sets  and  bases  for 
W  will  be  obtained  from  this  set  by  using  the  transformations  T  described  in  Theorem  2.26.  We 
let  (as  earlier)  V  denote  the  vertices  of  the  cube  [0..1/2]d.  We  continue  to  use  the  abbreviation 
V  :=  F\{0}. 

It  is  clear  that  the  space  Sl(rj)  is  generated  by  the  half-shifts  (tj(-  +  v))  y,  of  77,  and  this 
is,  indeed,  the  usual  starting  point  for  most  of  the  wavelet  constructions  now  in  the  literature. 
However,  we  will  show  below  that,  because  of  (3.5),  £*(77)  =  Sl(4>),  and  therefore  £l( 77)  =  «S(4>) 
with 

(3.6)  $  :=  (<£„  :=<£(•  +  t>))„€v 

The  generating  set  is  attractive  since  the  generator  <j>  of  S(<f>)  is  one  of  its  elements. 

Theorem  3.7.  If  £(<£)  C  (77)  and  supp  $  =  supp  rj,  then 

(3.8)  S'M  =S'(*)=  S(i) 
with  ■t  :=  (<K- +  ”)) 

Proof.  The  second  equality  is  clear.  As  for  the  first  equality,  let  g  :=  rj(-/2)  and  /  := 
<f>( -/2).  Since  <f>  6  Sl{r})  ,  we  have  /  €  S(g).  By  assumption,  supp g  =  supp/.  Therefore,  by 
Corollary  2.2,  £(/)  =  S(g).  Our  claim  then  follows  from  the  fact  that  Sx((f>)  and  Sl(ij)  are  the 
2-dilates  of  £(/)  and  S(g)  respectively.  6 

It  is  very  simple  to  find  elements  of  W.  If  /  €  S 1  (77),  then  since  P4,  is  the  orthogonal  projector 
onto  £(</>),  the  error  /  -  P^f  is  in  W.  If  we  choose  2d  -  1  such  functions  /  in  an  appropriate  way, 
we  obtain  a  basis  for  W.  Most  wavelet  constructions  begin  with  the  functions  /  =  g(-  +  v),  v  €  V . 
However,  there  are  too  many  of  these  functions  and  one  of  them  must  somehow  be  eliminated 
(destroying  symmetry).  Our  last  theorem  asserts  that  Sl(r])  is  also  generated  by  (<£(•  +  u))vgV. 
Starting  with  the  (<f>v  :=  <f>(-  +  v))v^v  gives  the  set 

W  :=  (wv  <f>v  —  P(p<f>v) V£yf 

It  is  easy  to  see  that  W  is  a  generating  set  for  W,  i.e.,  W  =  «S(W).  Indeed,  since  P<t>4>v  is  in  S(<f>), 
4>v  must  be  in  S{4> )  ©  S(wv).  Therefore,  {<£}  U  W  generates  Sx(t]).  It  follows  that  W  generates  W. 
From  (2.11),  we  obtain  a  simple  description  of  the  Fourier  transform  of  the  wj,: 

A. 

(3.9)  wv:=Tv-P^v  =  K-  »^':=n{0}. 

[ } 
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Theorem  3.10.  If  S(<f>)  C  S*( T7)  and  supp</>  =  supp  7?,  then  W  :=  Si  ©  <So  is  a  finitely  generated 
shift-invariant  space  and  W  :=  (wv)veV,  (defined  as  in  (3.9))  is  a  generating  set  for  W: 

(3.11)  W  =  <S(W). 

If  supp  4>  =  IRd,  then  W  provides  a  basis  for  W. 

Proof.  We  have  already  shown  that  W  is  a  generating  set  for  W .  If  supp  <f>  =  IR1*,  then  by 
Corollary  2.34,  the  set  $  provides  a  basis  for  <S($).  We  have  shown  that  :=  {< j> }  U  W  is  another 
generating  set  for  <S($).  Since  =  #$,  Theorem  2.26(ii)  asserts  that  also  provides  a  basis 
for  5($).  Hence,  det  <j($.)  is  nonzero  a.e.  From  the  orthogonality  between  W  and  we  find 

that  det  <?($«)  =  tfr  detG(VV).  Therefore,  det  G(W)  is  also  nonzero  a.e.,  thus  W  provides  a  basis 
for  W.  * 

The  set  W  is  our  first  set  of  generators  for  W.  We  shall  find  several  others  in  the  following 
sections.  The  idea  is  simple.  We  transform  W  to  a  new  set  of  generators  by  using  one  of  the 
matrices  T  whose  entries  are  2?r-periodic  functions.  The  intent  is  to  choose  T  in  such  a  way  that 
the  new  set  (TVV)V  of  generators  has  more  favorable  properties.  Our  next  result  illustrates  this 
procedure  and  provides  a  set  of  generators  which  are  compactly  supported  functions  whenever  4> 
is,  -  a  property  the  generators  in  W  lack. 

Theorem  3.12.  Assume  thatS(<f>)  C  Sl  (tj)  ,  supp  4>  =  suppq  =  IR  ,  and  that  (or equivalently 
4>)  is  bounded.  Then  the  2d  -  1  functions 

Wc:= 

provide  a  basis  for  the  wavelet  space  W .  If  <f>  has  compact  support,  then  the  functions  in  Wc  are 
also  of  compact  support. 

Proof.  Since  Wc  is  obtained  from  W  by  multiplying  by  the  27r-periodic  scalar  matrix 
[d>,  <t>]I,  which  is  assumed  here  to  be  bounded,  we  conclude  that  Wc  C  JDjCIR'1).  Furthermore, 

^  j 

suppc^  =  IR  implies,  by  Theorem  3.10,  that  W  provides  a  basis.  Hence,  by  Theorem  2.26,  Wc 
provides  a  basis  for  W  as  well. 

It  remains  to  show  that  the  functions  in  Wc  are  compactly  supported  whenever  <f>  is.  By  Lemma 
2.8,  [<l>v>4>]  is  a  trigonometric  polynomial.  Therefore,  the  inverse  transform  of  [4>v,<p\4>  is  a  finite 
linear  combination  of  the  shifts  of  <f>,  hence  is  compactly  supported  since  <f>  is.  The  same  argument 
shows  that  the  inverse  transform  of  [<f>,  <p](f>v  is  also  compactly  supported,  and  thus,  indeed,  the 
functions  in  Wc  are  compactly  supported.  4 

The  next  theorem  shows  that  W  always  has  a  set  of  generators  consisting  of  orthogonal 
wavelets. 

Theorem  3.13.  Let  S((f>)  C  S^tj)  and  supp  4>  =  supp  t)  =  IR'*.  Then: 

(i)  there  is  a  set  of  generators  9  for  W  which  provides  an  orthonormal  basis  for  W; 

(ii)  if  in  addition  <f>  has  compact  support,  then  there  is  a  subset  ^  =  (i Pv)v^v ,  of  compactly 
supported  functions  from  W  which  provides  a  basis  for  W  and  satisfies 

S(xl>v)  L  S(rl>u), 
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Proof.  (i):  By  Theorem  3.10,  W  is  a  basis  for  W  and  therefore  we  need  only  apply 
Theorem  2.28(i).  (ii):  By  Theorem  3.12,  the  functions  in  Wc  provide  a  basis  for  W  and  are  of 
compact  support  whenever  <f>  is.  Therefore,  we  need  only  apply  Theorem  2.28(ii).  6 

We  shall  next  discuss  conditions  under  which  the  half-shifts  w(-  +  v),  v  £  V',  of  a  function 
w  £  £2(1^)  provide  a  basis  for  W.  Clearly,  w  must  be  an  element  of  51(»7). 

Theorem  3.14.  Let  S((p)  C  S 1(rj)  and  supp0  =  fj  =  IRd.  Then: 

(i)  Ifw  £  Sl(rf),  then  the  functions  w(-+v),  v  £  V' ,  are  all  in  W  if  and  only  if[w,<j)]  is  2t -periodic. 

(ii)  If  w  is  a  generator  for  S1(tj)  and  [tu,  (f>\  is  2n -periodic,  then  the  functions  (w(-  +  v))  provide 
a  basis  for  W. 

Proof.  (i):  Since  w  €  L^IR,*),  the  function  is  in  £i(2'IT<<).  Proceeding  as  in  (2.6), 

(3.15)  f  w(x-  j)4>(x)dx  =  — 7x7^  /  e~j[w,4>\,  j  €  TLd  j2. 

Jnd  (*v)  Jind  \^r  J 2TT-1 

Now,  the  functions  w(-  +  v),  v  6  V' ,  are  all  in  W  if  and  only  if  the  inner  products  appearing  in 

(3.15)  are  zero  whenever  j  £  (2Zd/2)  \  7Zd,  i.e.,  if  and  only  if  [u5,<£]  has  period  27r. 

(ii):  From  the  facts  that  the  half-shifts  of  w  generate  Sl(rj),  and  supp?)  =  IRd,  we  easily 
conclude  that  supp  w  =  IR^,  and  therefore,  by  Corollary  2.34,  W.  :=  (^(•  +  1;))^^  provides  a  basis 
for  *S* (t?)-  Equivalently,  G(W.)  is  non-zero  a.e.  on  TT*.  It  follows  that  the  Gramian  of  any  subset 
of  W„,  and  in  particular  the  subset  W„\{u;},  is  non-zero  a.e.  as  well,  while,  by  (i),  this  latter  set 
lies  in  XV ,  since  we  also  assume  that  is  2x-periodic.  Thus,  we  have  found  2d  —  1  functions 

in  IV  (viz.,  the  functions  in  >V«\{io})  whose  Gramian  is  non-zero  a.e.,  while  Theorem  3.10  asserts 
that  W  contains  a  basis  of  cardinality  2d  -  1.  Thus,  Theorem  2.26(iii)  ensures  that  W„\{u>)  is 
basis  for  W .  4 

We  give  some  examples  of  functions  w  which  satisfy  the  assumptions  of  Theorem  3.14.  Since 
any  function  u;  £  W  is  in  £*(77),  we  must  have  w  =  rrj  for  some  47r-periodic  r,  and  so 

(3.16)  [ui,$J  =  [rrj,  4>]  =  r[fj,  4>\. 

We  would  like  the  function  in  (3.16)  to  be  of  period  2 n.  One  obvious  choice  is  to  take  r  =  1/(77,  <f>\ 
which  gives  the  function  wq  with  Fourier  transform 

(3.17)  wQ  :=  rj/[rj,ct)]  =  4>/[4>,4>\. 

The  half-shifts  u?o(-  +  v),  v  £  V' ,  will  generate  W  provided  Wq  €  ^(IR^)-  Note  the  intimate 
relation  between  the  present  wo  and  the  “dual  function”  p  which  was  defined  in  (2.15).  Indeed,  wq 

j  ^ 

is  orthogonal  to  each  4>( ■  -(-  j),  j  £  T2r  j 2\{0}  since  [«>o,  <i>\  —  1- 

In  general,  wq  will  not  be  in  Z-2 (IR-  )  because  of  the  division  by  =  <fr.  However,  we  can 

multiply  u?0  by  any  27r-periodic  function  and  obtain  the  Fourier  transform  of  other  candidates.  For 

example,  multiplying  by  ^(*  +  p)2  gives  the  function  w  with  Fourier  transform 

(3.18)  *  w:=4>  J]  ?(-  +  ^)2’ 
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This  function  w  also  has  the  advantage  of  being  of  compact  support  whenever  <j)  is.  Since  <f>  is  a 

generator  for  (77)  and  <f>2  >  0  a.e.  (because  supp<£  =  Hld),  the  function  w  of  (3.18)  is  a  generator 
for  Sl(ri)  (because  its  support  is  111“*).  Applying  (ii)  of  Theorem  3.14,  we  obtain  the  following 
Corollary. 

Corollary  3.19.  Let  S(4>)  C  Sl(q)  and  supp</>  =  supp rj  =  IR*.  If  the  function  w,  defined  by 
(3.18),  is  in  L2(IR‘*),  then  the  set  (ru(-  +  u)  :  c£  V'}  provides  a  basis  for  W.  If  4>  has  compact 
support,  then  so  does  w. 

We  note  that  the  above  w  is  in  I2  whenever  <f>  or  <t>  is  bounded.  For  example,  this  is  the  case 
whenever  <j>  has  compact  support. 

We  have  shown  so  far  that  it  is  easy  to  obtain  generating  sets  for  W  with  various  properties. 
They  can  be  chosen  to  provide  a  basis  or  an  orthogonal  basis  for  W.  They  can  be  chosen  to  be  the 
shifts  of  one  function  w  and  to  have  compact  support  if  <f>  has.  There  is  one  important  problem 
we  have  not  yet  discussed,  and  that  is  how  to  find  an  Z2(lR'*)-stable  basis  for  W  consisting  of 
compactly  supported  functions.  It  is  easy  to  see  that  the  generating  set  Wc  will  have  this  property 
if  the  half  shifts  of  are  Z,2(IRti)-stable  (this  assumption  is  not  realistic  in  the  stationary  case, 
but  can  be  satisfied  in  other  situations,  cf.  §8).  We  discuss  this  problem  in  §5  (in  the  univariate 
case)  and  in  §7  (for  the  multivariate  case).  But,  first,  we  examine  in  the  next  section  the  other  two 
conditions  in  (1.3),  i.e.,  (1.3)(ii,iii). 


4.  Multiresolution 

In  this  section,  we  analyze  conditions  (1.3)(ii)  and  (iii).  Our  setting  is  as  follows.  We  have 
for  each  k  £  7L  a  function  <f>k  £  Z/2(IRd)  and  the  space  Sk  Sk(4>k)  generated  by  its  2_fc-shifts. 
Alternatively,  Sk  is  the  2fc-dilate  of  the  principal  shift-invariant  space  «S(<At(2-fc-)).  We  noted  in 

(2.1)  that  the  functions  in  s  £  Sk  are  characterized  by  the  representation 

(4.1)  s  =  T<t>k,  t  of  period  2*+17t. 

We  first  study  condition  (1.3)(ii).  It  should  be  noted  that  we  have  completely  characterized  in 
[BDR]  density  properties  of  (arbitrary)  shift-invariant  subspaces  of  L2(IRd).  However,  the  present 
setting  is  so  simple  that  it  does  not  require  any  of  this  general  machinery. 

We  begin  with  the  following  lemma: 

Lemma  4.2.  Let  Sk,  k  £  7L,  be  a  nested  sequence.  Then  U Sk  is  a  closed  translation-invariant 
subspace  of  L2( fild). 

Proof.  Let  X  :=  L)Sk.  Then  X  is  certainly  closed.  Now,  let  f  £  X .  Since  Sk  C  Sk+1, 
f  £  Sk  for  all  k  sufficiently  large.  Since  Sk  is  2-fc-shift-invariant,  ft  /(•  +  t)  is  in  X  for  any 
t  =  2 ~kj,  j  £  7Zd,  which  means  that  ft  is  in  X  for  all  dyadic  t  =  2 ~kj,  j  £  7Ld,  k  £  TL.  Since 
translation  is  a  continuous  operation  in  L2(IRd),  we  obtain  that  ft  is  in  X  for  all  t  £  IRd.  Moreover, 
if  g  £  X,  and  f  £  X,  then  ||sf(  -  /tj|  =  ||<7  -  f\\.  Approximating  g  by  functions  f  £  X  shows  that 
9t  €  A.  * 
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It  is  well  known  that  a  closed  translation-invariant  subspace  X  of  £2(IRd)  is  characterized  by 
its  Fourier  transforms.  Precisely,  X  =  £2(0)  for  some  measurable  set  f 1  (called  the  spectrum  of 
X). 

Theorem  4.3.  Let  (Sk  :=  be  a  nested  sequence.  Then  USk  =  L2(JR-d)  if  and  only  if 

fl0  :=  Usupp^>fc  =  JR*  ( modulo  a  null-set). 

Proof.  Let  X  :=  U Sk.  From  the  above  remarks  on  translation-invariant  spaces,  X  = 
£2(Q)  for  some  measurable  set  Cl  C  IRd.  We  have  X  =  £2(IRd)  if  and  only  if  0  =  IR"1  modulo  a 
null-set.  We  verify  that  Cl  =  CIq  modulo  a  null-set  which  will  complete  the  proof.  Since  each  <j>k 
is  in  X,  we  must  have  supp^t  C  Cl  modulo  a  null-set,  and  so,  CIq  C  Cl  modulo  a  null-set.  Now 
suppose  that  Cl\Cl0  contains  a  set  Cli  of  positive  measure.  From  (4.1),  each  element  in  Sk,  k  £  ZZ, 
has  Fourier  transform  which  vanishes  on  fii.  Hence,  each  element  in  U Sk  has  Fourier  transform 
which  vanishes  on  flj.  Taking  the  closure,  we  see  that  each  element  in  X  has  Fourier  transform 
which  vanishes  on  Cl j.  This  is  absurd  since  X  contains  L2(Cli).  4 

The  role  of  (1.3)(ii)  in  multiresolution  analysis  is  to  guarantee  that  lim/t-.oo  Pkf  =  /  for  each 

/  €  £a(ffi*). 

Corollary  4.4.  Let  ( Sk  :=  Sk(<j>k)) ke7l  be  a  nested  sequence,  and  UfcezsuPP  <t>k  =  IRd.  Then, 
the  orthogonal  projectors  Pk  from  £2(IR<3t)  onto  Sk  satisfy  limfc_oo  Pkf  =  /  for  all  f  £  L2(lR.d). 

Proof.  Since  Sk  C  Sk+1 ,  Theorem  4.3  says  that  ||/  -  Pkf\\  =  dist(/,5A:)  — ►  0.  A 

We  next  consider  in  more  detail  the  stationary  case  Sk  =  «Sfc(0(2fc-)),  k  £  ZZ,  i.e.,  the  case 
when  Sk  is  the  2fc-dilate  of  «S°  =  S(4>),  which  is  the  usual  situation  treated  in  multiresolution.  The 
following  is  a  very  simple  sufficient  condition  for  (1.3)(ii),  in  the  event  of  a  stationary  multiresolu¬ 
tion. 

Theorem  4.5.  Let  </>  be  an  £2(IRd)-function,  and  for  each  k  £  TL,  let  Sk  be  the  2k-dilate  of 
Assume  that  (Sk) .  is  nested.  Then,  ( Sk)k  satisfies  (1.3)(ii)  if  4>  is  non-zero  a.e.  in  some 
neighborhood  of  the  origin. 

I  ^  ^ 

Proof.  Here  (f>k  =  <£(2  •),  and  therefore  <pk  =  Ck<t>(-/ 2").  Thus  if  <f>  is  non-zero  a.e.  on  Cl, 
then  <t>k  is  non-zero  a.e.  on  2  Cl.  Now,  if  Cl  is  some  neighborhood  of  the  origin,  we  obtain  that 
U*supp<^t  =  Dld,  since  Ufc2fcQ  =  Dld.  By  Theorem  4.3,  (1.3)(ii)  holds.  4 

Of  course,  (1.3)(ii)  can  also  hold  when  <f>  vanishes  at  every  neighborhood  of  the  origin  on  a  set 
of  positive  measure.  For  example,  this  is  the  case  if  d  =  2  and  supp  ^  =  {16  IR2  :  A  <  |*i|}. 

Special  cases  of  Theorem  4.5  have  been  established  by  other  authors.  For  example,  in  the 
univariate  case  and  under  certain  restrictions  on  the  smoothness  and  decay  of  <f>,  Mallat  [Ma] 
showed  that  whenever  <p  has  orthonormal  shifts,  assumption  (1.3)(i)  implies  that  \jSk  =  L2{lRd). 
Recently,  this  was  generalized  to  the  multivariate  case  by  Jia  and  Micchelli  [JM]  who  replaced 
orthonormality  by  £2(IRrf)-stability  and  replaced  Mallat’s  other  conditions  by  the  requirements 
that  <f>°  £  L2(TLd),  and  that  4>  satisfy  the  refinement  equation  (1.6)  for  a  sequence  a  from  Ci(7Zd). 
In  both  of  these  examples,  the  conditions  used  imply  that  <j>  is  continuous  and  <f>( 0)  ^  0,  hence 
these  results  indeed  follow  from  Theorem  4.5.  Note  that,  by  the  same  token,  Theorem  4.5  certainly 
applies  whenever  <j>  £  L\(\fld)  H  £2(IR<1)  and  <£(0)  ^  0. 
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We  return  to  the  general  case  of  the  spaces  Sk  =  Sk(<j>k )  introduced  at  the  beginning  of  this 
section,  in  order  to  discuss  condition  (1.3)(iii).  We  shall  need  the  notion  of  Lebesgue  points:  recall 
that  if  /  is  locally  in  L\{ IIId),  a  point  x  £  IRd  is  said  to  be  a  Lebesgue  point  of  /  if 

in.mnI7)i  /  /(“)  du  =  /(*)> 

IQI-o  IQ]  JQ 

with  the  limit  taken  over  cubes  Q  which  contain  x.  For  each  locally  integrable  /,  almost  every 
point  is  a  Lebesgue  point  (see  e.g.  [BS;  p.  121]).  We  shall  need  the  following  simple  lemma  (which 
certainly  is  known). 

Lemma  4.6.  If  Cl  is  a  measurable  subset  oflRd  and  a  ^  0  is  a  fixed  real  constant  such  that  for  each 
dyadic  t  £  IRd,  we  have  f l  +  at  =  Cl  modulo  a  null-set,  then  Cl  =  IRd  or  Cl  =  0,  modulo  a  null-set. 
Moreover,  if  f  is  any  measurable  function  on  IRd  which,  for  each  dyadic  t,  satisfies  /(•  -f  at)  =  / 
a.e.  ,  then  f  =  const  a.e. 

Proof.  By  replacing  /  (respectively,  Cl)  by  /(a-)  (respectively,  a'1!!),  we  can  assume  that 
a  =  1.  Let  a  £  IRd  be  a  Lebesgue  point  of  \  '•=  Xn-  Then  with  Q$  :=  [-6/2 .  .6/2\d,  we  have 

(4.7)  Yim  6~d  f  x{x  +  a)dx  =  x(a). 

5-°  JQi 

Now,  for  any  dyadic  number  t,  we  have  x(*  +  <)  =  x(x)>  a-e-  *n  x-  Hence,  for  any  set  Q  of  finite 
measure, 

(4.8)  f  x(z  +  a  +  t)  dx  =  /  x(*  +  aMx- 

Jq  Jq 

If  y  €  Htrf  is  any  other  Lebesgue  point  of  x,  then  using  the  density  of  the  dyadic  points,  we  can,  for 
each  6  >  0,  find  a  ts  such  that  y  €  a  +  ts  +  Qs .  Using  (4.8),  we  find 

X(y)  =  lim  6~d  [  x(x  +  a  +  tg) dx  =  lim  6~d  f  x(x  +  a) dx  =  x(a)- 
«-°  Jq,  6-+o  JQf 

Hence,  x  is  constant  a.e.  and  our  result  follows. 

If  the  function  /  is  as  described  in  the  lemma’s  statement,  then  for  each  y  €  IR,  the  set 
Cl  :=  {x  :  /(*)  <  y}  satisfies  Cl  +  t  =  Cl,  modulo  a  null-set,  for  each  dyadic  t.  Hence  11  =  lRd 
or  Cl  =  0,  modulo  a  null-set,  and  it  follows  that  /  is  a  constant  a.e.  If  /  is  complex-valued,  this 
argument  can  be  applied  to  both  its  real  and  its  imaginary  parts.  4 

Theorem  4.9.  Given  the  sequence  Sk  :=  Sk(<f>k),  k  £  TL,  set  Y  :=  C\kSk.  Then  Y  is  a  linear 
subspace  of  L 2(nid)  of  dimension  <  1. 

Note  that  the  nestedness  assumption  (1.3)(i)  is  not  made  here.  Further,  as  the  proof  below 
makes  clear,  the  result  remains  valid  even  when  {«Sfc}fc  is  replaced  by  a  subsequence  {«Sfc»  }j,  provided 
that  lim j _ oo  kj  =  -oo. 
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Proof.  Assuming  Y  ^  {0},  we  will  show  that  dimF  =  1.  Let  f,g  be  two  functions  in  Y. 
We  define  the  following  function 


F  :  IR^  -+  C2  :  i  h 


(/(j).g(j)) 

ll(/(*).7(*))lf 


f(x)  =  g(x)  =  0, 
otherwise. 


We  shall  prove  that  F  is  constant  on  its  support.  Let  if  be  a  measurable  subset  of  C2\{0},  and 
assume  that  A  :=  F~l(I\)  has  positive  measure.  Let  D  be  the  set  of  all  points  of  the  form  2kjir, 
k  €  2Z,  j  €  2Zd.  Let  x  £  A,  and  t  =  2k+1jir  £  D.  In  view  of  (4.1),  there  are  r  and  r)  of  period 
2fc+1ir  such  that 


(4.10)  f  =  r$ky  g  =  T)4>ky  a.e. 

Since  F(x)  ^  0  (since  F(x)  £  K  and  K  does  not  contain  0),  4>k{x)  ^  0,  and  we  have 

(/(*  +  t),g(x  +  t ))  =  <i>k(x  +  t)(r(z  +  0*  »?(*  +  0)  =  $k(x  +  t)(r{x),  tj(x))  =  -4^ -~(f(x),g{x)). 

<t>k(z) 

This  implies  that  either  F(x  +  t)  =  0  or  that  F(x)  =  F(x+t).  We  conclude  that^A+Z?)  C  A'u{0}. 
Since  A  is  assumed  to  be  of  positive  measure,  so  must  be  A  4-  D,  and  hence  Lemma  4.6  (when 
applied  to  U  :=  A  +  D)  implies  that  A  +  D  =  Ilt^  (up  to  a  null-set).  Thus  (again  up  to  null  sets)  F 
assumes  its  non-zero  values  in  K ,  and  since  K  can  be  made  arbitrary  small,  this  can  happen  only 

>"V. 

if  F  assumes  essentially  only  one  non-zero  value.  Equivalently,  the  functions  /  and  g  are  linearly 
dependent,  which  is  what  we  wanted  to  prove.  4 

It  can  indeed  happen  that  Y  has  dimension  1.  For  example  if  <f>k  =  <j)  for  each  k  (not  to  be 
confused  with  the  stationary  case:  <pk  =  <t>(2k-),  all  k),  then  <f>  is  obviously  in  Y .  Other,  less  trivial, 
examples  are  also  possible  (see  §6).  In  passing,  we  note  the  following  immediate  consequence  of 
Corollary  2.2  (and  its  scaled  versions): 

Proposition  4.11.  Let  f  £  C\kzTZ.Sk(<t>k)-  Then  f  generates  all  the  spaces  Sk  if  and  only  if 
supp  f  =  supp 4>kl  all  k.  In  particular,  the  spaces  (Sk(<t>k))k  are  generated  all  by  a  single  function 
only  if  supp  <f>k  =  supp^/  for  all  k,k'  £  TL. 

We  also  note  that,  for  any  nested  sequence  Xk  C  of  closed  subspaces  of  a  Hilbert  space, 
the  corresponding  orthogonal  projectors  Pk  :=  Pxk  converge  strongly  to  the  orthogonal  projector 
P- oo  onto  A_oo  :=  flfeX*  as  k  — »  —00  (hence  converge  strongly  to  the  orthogonal  projector  Poo 
onto  Xoo  :=  UfcAfc  as  k  — ►  00).  Therefore,  in  particular: 

Theorem  4.12.  Let  Sk  :=  Sk(4>k),  k  £  TL,  be  a  nested  sequence,  and  let  Y  =  nk£’zSk  be  the 
(one-  or  zero-dimensional)  space  of  Theorem  4.9.  Then  limjt__oo  Pkf  =  Pyf,  for  all  f  £  L2(IR<i)- 

Proof.  Here,  for  completeness,  is  a  proof  which  only  uses  the  fact  that  Xk  :=  Sk  is  a 
nested  sequence  of  closed  linear  subspaces  of  a  Hilbert  space,  with  Y  :=  C\kXk. 

To  show  that  linu__oo  Pkf  =  Pyf,  it  is  sufficient  to  show  that  Pkf  g  weakly  implies  that 
g  =  Pyf  and  g  =  limfc__oo  Pkf.  For,  it  implies  that  Pyf  is  the  only  limit  point  of  (Pkf)k,  and, 
further,  implies  that  (Pkf)k  has  limit  points  since,  being  bounded,  it  has  weak  limit  points. 
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So,  let  g  be  the  weak  limit  of  ( Pkf)k •  Since  every  Xj  is  closed  and  convex,  hence  weakly 
closed,  and  contains  every  P*/  with  k  <  j,  it  must  contain  g;  therefore,  g  £Y.  On  the  other  hand, 
x*  :=  f  —  Pkf  is  perpendicular  to  Xk,  hence  to  Y,  therefore,  so  is  the  weak  limit,  x  :=  f  —  g.  In 
short,  g  =  Pyf-  Since  Y  is  the  intersection  of  the  nested  sequence  (Xk) k,  it  follows  that 

lim  11**11  =  sup dist(/,  Xk )  <  dist(/,  Y)  =  ||x||. 

fc— »  — oo  k 

Because  of  the  weak  convergence  we  also  have  ( x*,x )  — *  ( x,x ),  and  therefore 

0  <  || Pkf  -  0||2  =  ||*  -  **||2  =  ||x||2  -  2  Re(x,  xk)  +  ||**||2  —  -||x||2  +  lim  ||x*||2  <  0 

k 

as  k  runs  to  -oo.  In  other  words,  lim*  P*/  =  g  -  Pyf-  4 

Since  Y  C  Sk,  it  is  orthogonal  to  each  of  the  wavelet  spaces  Wk  :=  Sk+1  ©  Sk,  k  £  TL.  There¬ 
fore,  applying  Corollary  4.4  and  Theorem  4.12  ,  we  obtain  the  following  orthogonal  decomposition 
of  L2(lRd). 

Corollary  4.13.  Let  Sk  :=  Sk(<f>k ),  k  £  7L,  be  a.  nested  sequence,  and  let  :=  Usupp^*  =  IRd 
(modulo  a  null-set).  IfY  is  the  (one-  or  zero-dimensional)  subspace  of  Theorem  4.9,  then 

l2(  iRrf)  =  y  ©  ©w*. 

kezi 

The  significance  of  the  last  corollary  is  the  following.  Let  ifry  be  any  nontrivial  element  of  Y 
with  ||0k||  =  1.  If  for  each  k  £  7L,  the  set  provides  an  orthonormal  basis  for  the  wavelet  space 
YVk,  then  the  totality  of  functions  rjjy  and  0(-  -  j 2~k),  j  £  7Ld,  rp  £  k  £  TL ,  is  an  orthonormal 
basis  for  LtCHT*).  Thus,  even  when  Y  is  nontrivial,  multiresolution  produces  a  basis  for  L2(lR<i). 
Similarly,  we  obtain  an  L2-stable  basis  whenever  the  ’P*  provide  an  Z,2-stable  basis  for  Wk  whose 
stability  constants  are  independent  of  k  £  TL. 

In  the  stationary  case,  i.e.,  the  case  when  =  ^>o(2fc-),  the  following  corollary  shows  that  Y 
is  necessarily  trivial. 

Corollary  4.14.  For  4>  £  /^(nr1),  define  Sk  :=  Sk {<f>(2k ■)) ,  k  £  TL.  Then  f)Ske2z  =  {0}. 

Note  that,  as  in  Theorem  4.9,  the  nestedness  condition  (1.3)(i)  is  not  required,  hence  is  not 
assumed. 

Proof.  We  suppose  that  /  is  a  non-trivial  function  in  T2(IR<i)  which  is  in  each  of  the  spaces 
Sk  and  derive  a  contradiction.  By  the  assumptions  here,  each  Sk  is  the  2-dilate  of  Nfc_1,  and  hence 
n*«Sfc  is  invariant  under  dilation  by  2.  On  the  other  hand,  by  Theorem  4.9,  this  space  is  at  most 
one-dimensional.  Therefore,  if  f  £  then  there  exists  some  A  such  that 

(4.15)  /(2  )  =  A/,  a.e.  on  IRd. 

It  is  now  easy  to  show  that  this  is  impossible  for  /  6  L2(IR<1)\{0}.  Indeed,  for  each  C  >  0,  the  sets 
Fk  {x  :  2k  <  |x|  <  2fc+1  and  |/(x)|  >  C|A|fc}  satisfy 

Fk  —  2Fk-\  and  meas(F*)  =  2dmeas(Ffc_i),  for  all  k. 
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If  /  is  not  the  zero  function,  then  for  some  C  >  0,  mea s(F0)  ^  0.  From  (4.15),  |/(x)|  >  C\k  for 
x  €  2 kF0.  Therefore, 

/  |/|!>C!meas(/„)^(2J|A|J)fc, 

JlRd  *€2Z 

which  shows  that  /  is  not  in  L 2(1R<1)  because  the  series  diverges.  4 

The  importance  of  Corollary  4.14  is  that,  in  the  stationary  case,  it  is  not  necessary  to  assume 
property  (1.3)(iii).  Moreover,  in  the  case  that  <f>  has  compact  support,  condition  (1.3)  (ii)  is  already 
implied  by  (1.3)(i).  We  have  therefore  the  following  corollary. 

Corollary  4.16.  If,  for  the  compactly  supported  function  <f>  in  Z,2(IRd),  the  sequence  Sk  := 
Sk(<f>( 2k-)),  keTL,  is  nested,  i.e.,  satisfies  (1.3)(i),  then  the  conditions  (1.3)(ii,iii)  are  automat¬ 
ically  satisfied  and  we  have  the  orthogonal  decomposition 

l2(  nt*)  =  ©  wk 

k<m 

with  W  :=  Sl  ©  S°  the  wavelet  space  and  Wk  its  2k-dilate,  k  6  ZZ. 


5.  Univariate  wavelets  and  prewavelets 

After  showing  in  the  last  section  that  conditions  (ii)  and  (iii)  of  (1.3)  hold  in  quite  a  general  set¬ 
ting,  we  now  turn  our  attention  back  to  wavelet  constructions.  We  start  with  a  separate  discussion 
of  the  univariate  case,  since  this  case  is  significantly  simpler  than  its  multivariate  counterpart. 

As  in  §3,  we  are  only  interested  in  studying  one  of  the  wavelet  spaces,  namely,  W  :=  S1  Q  S°. 
The  other  wavelet  spaces,  Wk  :=  «S*+I  0  Sk ,  are  obtained  by  identical  methods,  and  furthermore, 
in  the  stationary  case  each  of  the  wavelet  spaces  is  obtained  from  W  by  dilation. 

We  work  in  the  same  setting  as  in  §3:  We  assume  that  <f>,  q  6  ^(IR-)  satisfy 

(5.1)  S{4>)CS\r,) 
and 

(5.2)  supp  <f>  =  supp  rj  =  IR. 

As  before,  we  remind  the  reader  that  this  last  assumption  is  always  satisfied  when  4>  and  77  are 
compactly  supported. 

We  have  seen  in  §3  that  W  is  a  principal  shift-invariant  space  and  there  is  always  a  generator 
w,  for  W  whose  shifts  form  an  orthonormal  system.  However,  in  general  we  can  say  nothing  about 
the  support  of  in.,  or  the  decay  of  u?.(x)  as  ]x|  — *  00.  In  this  section,  we  want  to  go  further  and 
find  other  generators  for  W  with  favorable  decay  properties.  In  particular,  in  the  case  rj  =  4>{ 2-) 
usually  studied  in  wavelet  constructions,  we  shall  recover  various  generators  for  W  given  by  other 
authors. 
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Our  starting  point  is  the  function  w  of  Theorem  3.10,  i.e.,  the  function  whose  Fourier  transform 


is 


(5.3) 


w  :=  e1/2£- 


- - <p. 

(<£,  $ 


We  know  from  Theorem  3.10  that  w  provides  a  basis  for  W. 

We  wish  to  express  the  generators  that  follow  in  terms  of  rj.  For  this  we  shall  use  the  refinement 
relation  (1.6): 


(5.4) 


4>~  At) 


with  A  a  4T-periodic  function. 

Theorem  5.5.  Let 

(5.6)  :=  2e_i /2A(-  +  2t r)^(-  +  2 k)2tj  =  2e_i/2[f],<t>](-  +  2k)t}. 

If  p  €  X-2(IR.),  then  its  inverse  transform  p  is  a  generator  for  the  wavelet  space.  Moreover,  p  has 
orthonormal  ( respectively ,  stable)  shifts  if  the  shifts  of  p  and  the  half-shifts  of  p  are  orthonormal 
(respectively,  stable). 

Proof.  We  already  know  that  the  function  w  of  (5.3)  is  a  generator  for  W.  We  will  show 

/'■v 

that  w/p  is  2ff-periodic.  Since  (5.2)  implies  that  supp  p  =  HI,  this  will  prove,  by  (1.11),  that 
w  £  S(p)  and  hence  p  generates  W. 

Since  ei/2(-+27r)  =  -ei/2,  and  ex/2  is  47r-pcriodic,  we  see  that  [e1/20,$]  =  eX/2(p2 -<j>( -+27r)2), 
while  [p,p]  =  p2  +  p(-  +  2k) .  Substituting  this  into  (5.3),  we  obtain  that 


~  2<f>(-  +  2?r)2  ~ 

w  =  -  - e1/2tf. 

9* 


Since  p  —  At},  and  A  is  47r-periodic,  p  =  M|^-  Therefore,  we  see  that 

r  :=  w/p  =  ( p)~2A  /!(•  +  27r)ei. 

Since  A  is  4^-periodic  and  p  and  ex  are  27r-periodic,  we  conclude  that,  indeed,  the  ratio  w/p  is 
2x-periodic  and  hence  p  generates  W. 

To  prove  the  rest  of  the  theorem,  we  first  compute  p  as  follows: 

P2  =  4  {|/1(-  -f  27r)|277(-  +  27r)47 f  +  |/1|V^(-  +  2k)2  J 

(5.7)  =  4’H-  +  2k)2  r/2  +  2tt)|27?(-  +  2k)2  +  |4|2rp} 

=  4r)(-  +  2k)2t}2  |$(.  +  2k)2  +  p2  j  =  4t ;(•  +  2k)2t)2P2. 

If  p  has  orthonormal  shifts  and  r)  has  orthonormal  half-shifts,  then  <f>2  =  1  a.e.,  and  rj2  =  1/2  a.e. 
We  conclude  that  p2  =  1  a.e.,  and  hence  p  has  orthonormal  shifts.  Similarly,  if  p  has  T2(Dl)-stable 
shifts  and  r)  has  T2(IR)-stable  half-shifts,  then  the  functions  p  and  rj  and  their  reciprocals  are 
bounded.  It  follows  that  p  has  the  same  property  and  hence  the  shifts  of  p  are  T2(Dl)-stable.  4 
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Remark  5.8.  It  also  follows  from  (5.7)  that  the  stability  constants  C\{ip),C2(ip)  >  0  for  ip,  /.e., 
the  positive  constants  in  the  inequality  C\(ip)  <  ip  <  Cjiip)  a.e.,  can  be  chosen  as 

Cj(iP)  =  2Cj(<P)Cj(q)2,  j  =  1,2, 

where  Cj(<£)  and  Cj(q),  j  =  1,2,  are  the  stability  constants  associated  with  <p  and  q  respectively. 

Note  that  when  q  =  \/20(2-),  the  orthogonality  assumption  or  the  stability  assumption  on  the 
half-shifts  of  q  is  equivalent  to  the  corresponding  assumption  on  the  shifts  of  <p.  Further,  in  the 
orthogonal  case,  q( •  +  2tt )2  =  1/2,  hence  the  formula  (5.6)  is  reduced  to 

(5.9)  xp  =  e_l/2A(-  +  2ir)q, 

which  gives  the  usual  wavelet  obtained  by  multiresolution.  Note  also  that  the  theorem  incidentally 
proves  that  ip  €  £.2(IR.),  hence  ip  is  in  ^(Di)  whenever  q  has  L2(IR)-stable  half-shifts. 

Mallat  has  proved  the  orthonormal  part  of  the  above  theorem  (for  q  :=  y/2<p(2-))  without 
the  assumption  (5.2),  but  with  additional  hypotheses  on  the  decay  and  smoothness  of  <p.  Several 
authors  have  used  Mallat ’s  approach  to  construct  orthonormal  wavelets,  including  Daubechies  [D] 
in  her  celebrated  construction  of  wavelets  w  of  compact  support  and  arbitrary  high  orders  of  dif¬ 
ferentiability.  However,  the  difficult  part  of  the  Daubechies  construction  is  to  show  the  existence  of 
compactly  supported  functions  <p  which  satisfy  (5.4),  have  arbitrarily  high  orders  of  differentiability, 
and  have  shifts  which  are  orthonormal. 

As  an  example,  if  <p  is  the  B-spline  <p  =  Af(|0, . . .,  r)  of  order  r  with  knots  at  0,  ...,r,  then 
S(d>)  is  the  space  of  all  cardinal  splines  of  order  r  which  are  in  Z,2(Hf).  The  function  ip  is  then  the 
spline  wavelet  of  Battle  and  Lemarie  (see  [B]). 

The  prewavelet  part  of  Theorem  5.5  has  been  proved  by  Micchelli  in  [Mi],  but  under  different 
hypotheses.  He  docs  not  assume  (5.2),  but  assumes  that  <p  satisfies  the  refinement  equation  (1.6) 
with  coefficients  a  £  Similar  ideas  have  been  employed  by  Chui  and  Wang  [CW],  [CWl]. 

In  particular,  when  <p  is  the  cardinal  B-spline,  the  prewavelet  ip  of  Theorem  5.5  is  their  compactly- 
supported  spline  wavelet  (except  for  an  integer  shift). 

The  remainder  of  this  section  will  be  devoted  to  the  important  case  when  the  functions  <p  and  q 
are  compactly  supported.  We  shall  be  interested  in  finding  functions  w  from  W  which  have  minimal 
support.  In  the  context  of  compact  support,  one  encounters  the  notion  of  linear  independence:  We 
say  that  the  shifts  of  the  compactly  supported  <p  are  linearly  independent  if,  for  each  sequence 
c,  the  sum  <P( •  -  j)c[j)  is  identically  zero  if  and  only  if  c(j)  =  0  for  all  j  £  TL  (note  that  for 

each  x  €  IR,  the  series  has  only  a  finite  number  of  nonzero  terms  and  hence  converges  pointwise). 
We  remark  that  linear  independence  of  the  shifts  of  <p  implies  that  these  shifts  are  Z,2(IR)-stable, 
(cf.  [JM]). 

In  what  follows,  we  denote  by  diam  D  the  length  of  the  smallest  interval  containing  the  subset 
fi  of  HI.  With  the  aid  of  [R2],  the  following  result  on  linearly  independent  generators  was  proved 
in  [BDR1] . 

Result  5.10.  Let  S  be  a  univariate  principal  shift-invariant  space  which  is  generated  by  a  com¬ 
pactly  supported  function.  Then,  there  exists  a  compactly  supported  <p  £  $  that  satisfies  all  of  the 
following  conditions: 
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(a)  The  shifts  of  4>  are  linearly  independent. 

(b)  Every  compactly  supported  f  €  S  can  be  written  as  a  finite  linear  combination  of  the  shifts 
of  <j>. 

(c)  diam  supp  <£  <  diamsupp  /,  for  every  f  £S. 

Furthermore,  up  to  a  shift  and  a  scalar  multiplication,  4>  is  characterized  by  any  of  these  three 
properties. 

Corollary  5.11.  If  <f>  has  compact  support,  then  the  wavelet  space  W  has  a  compactly  supported 
generator  whose  shifts  are  linearly  independent.  This  generator  enjoys  all  the  properties  of  <f>  in 
Result  5.10. 

Proof.  By  the  case  d  =  1  in  Theorem  3.12,  W  is  principal  and  has  a  compactly  supported 
generator.  It  is  therefore  enough  to  apply  Result  5.10.  4 

In  view  of  the  attractive  properties  of  a  linearly  independent  generator,  it  is  desirable  to  find 
a  constructive  method  to  find  the  linearly  independent  generator  of  W .  For  this,  we  shall  assume 
(without  loss  of  generality  in  view  of  Result  5.10)  that  the  generator  77  for  5!( 77)  has  linearly 
linearly  independent  half-shifts.  In  view  of  Result  5.10,  any  compactly  supported  function  in  Sl(rj) 
has  Fourier  transform  rrj  with  r  a  4jr-periodic  trigonometric  polynomial.  We  are  interested  in  the 
properties  of  r  that  characterize  linear  independence  of  the  shifts  of  (r^)v. 

If  r  is  a  non-trivial  47r-periodic  trigonometric  polynomial,  then  r  =  em/2  Sj=o  w*th 

o(0)a(rc)  ^  0  and  m  €  21.  We  call  n  the  modified  degree  of  r  and  write  mdegr  :=  n.  From 
this  it  easily  follows  that  if  /  and  g  are  compactly  supported  and  g  =  tJ  for  some  4rr-periodic 
trigonometric  polynomial,  then 

(5.12)  diam  supp  g  -  diam  supp  /  +  (mdeg  r)/2. 

If  r  and  C,  are  two  47r-periodic  trigonometric  polynomials,  we  say  that  £  divides  r  if  t/C,  is  also 
a  4tt- periodic  trigonometric  polynomial.  With  this,  we  have  the  following  characterization  of  the 
linearly  independent  generators  of  W. 

Proposition  5.13.  Assume  77  is  compactly  supported  and  has  linearly  independent  half-shifts. 
Let  w  be  any  compactly  supported  generator  of  the  wavelet  space  W .  Then  w  is  the  linearly 
independent  generator  of  the  wavelet  space  (and  thus  enjoys  all  the  properties  of  the  <p  in  Result 
5.10)  if  and  only  if  the  Air-periodic  trigonometric  polynomial  r  in  the  representation  w  =  rrj  is  not 
divisible  by  a  non-constant  2ir-periodic  trigonometric  polynomial. 

Proof.  By  Result  5.10,  every  compactly  supported  w  £  W  C  Sx{g)  can  be  expressed  in 
the  form  u;  =  rrj  for  some  47r-periodic  trigonometric  polynomial  r.  If  r  —  A£,  where  A  is  a  2 it-  and 
(  is  47r-periodic  trigonometric  polynomial,  then,  by  Corollary  2.2,  w'  :=  (C^)v  is  also  a  compactly 
supported  generator  for  W. 

By  (5.12), 

diam  supp  w  =  diam  supp  77  +  (mdeg  r)/2,  diam  supp  w'  =  diam  supp  77  +  (mdeg£)/2, 

hence  diam  supp  <  diam  supp  u;,  with  equality  if  and  only  if  mdeg(  =  mdegr,  i.e.,  if  and  only 
if  mdeg  A  =  0.  Our  claim  follows  then  from  the  fact  (cf.  Result  5.10)  that  the  linearly  independent 
generator  is  characterized  by  the  minimality  of  its  support.  4 
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In  view  of  the  last  result,  the  search  for  the  linearly  independent  generator  of  W  can  be  carried 
out  as  follows:  assuming  77  has  linearly  independent  half-shifts,  we  find  some  particular  compactly 
supported  w  £  W,  and  write  w  =  rfj.  Then  r  is  necessarily  a  trigonometric  polynomial.  Factoring 
r  =  A£  where  A  is  a  27r-periodic  factor  of  maximal  degree,  ip,  defined  by  ip,  =  (rj  is  the  linearly 
independent  generator  of  W. 

Corollary  5.14.  Assume  that  <p  and  77  are  compactly  supported  and  the  half-shifts  of  77  are  linearly 
independent,  and  that  <p  is  compactly  supported,  and  <p  Arj.  Then  the  linearly  independent 
generator  ip,  for  the  wavelet  space  W  =  Sl(rj)  Q  S(<p)  is  given  by 

&  :=  c*?, 

where  (  :=  r/A,  and  A  is  a  2tt -periodic  trigonometric  polynomial  of  maximal  degree  that  divides 

t  :=  e_i/2^(-  +  2x)t}(-  +  2?r)2  =  e_i/2[»?,$|(-  +  2*). 

Proof.  From  Theorem  5.5,  we  know  that  ip  :=  (tt?)v  generates  W .  Thus  the  claim  will  fol¬ 
low  from  the  argument  preceding  this  corollary,  as  soon  as  we  show  that  r  indeed  is  a  trigonometric 
polynomial. 

The  function  rj2  is  a  trigonometric  polynomial  by  the  analogue  of  Lemma  2.8,  since  it  is  the  4ir- 
periodization  of  |^|2  for  the  compactly  supported  tj.  The  mask  A  is  also  a  trigonometric  polynomial 
by  Result  5.10,  since  t]  and  <p  are  compactly  supported  and  77  is  a  linearly  independent  generator 
oiSx(rj).  * 

Result  5.10  tells  us  that  the  search  for  a  linearly  independent  generator  is,  necessarily,  the  same 
as  the  search  for  a  minimally  supported  generator  in  the  sense  that  we  are  minimizing  diam  supp  w 
among  all  generators  w.  Chui  and  Wang,  [CW],  considered  a  slightly  different  notion  of  minimality: 
they  were  interested  in  finding  a  generator  w  for  W  which  can  be  expressed  in  the  form  w  =  rrj,  with 
r  a  trigonometric  polynomial  of  minimal  degree  (they  assume  that  the  refinement  mask  A  =  <p/rj 
is  a  polynomial,  to  guarantee  the  existence  of  such  r).  Thus,  while  we  minimize  diam  supp  w  over 
all  possible  generators  w,  Chui  and  Wang  minimize  diam  supp  w  only  over  those  w  which  can  be 
written  as  a  finite  linear  combination  of  the  half-shifts  of  77.  However,  because  of  Result  5.10,  the  two 
notions  coincide  if  we  assume  (as  we  do)  that  the  half-shifts  of  77  are  linearly  independent.  In  any 
event,  with  straightforward  modifications,  the  arguments  used  in  Proposition  5.13  and  Corollary 
5.14  can  be  applied  to  show  that  the  same  characterization  holds  for  the  “minimal  w"  in  the  [CW] 
sense. 

Chui  and  Wang  stated  their  results  in  terms  of  the  symmetric  zeros  of  the  trigonometric 
polynomials  involved.  Let  us  pause  for  a  moment  to  see  how  symmetric  zeros  enter  into  the 
characterizations  provided  above.  If  r  is  a  47r-periodic  trigonometric  polynomial,  then,  up  to  some 
exponential  factor,  we  can  write  r  =  p(e j/2)  for  some  algebraic  polynomial  p  with  degp  =  mdegr. 
But,  for  any  algebraic  polynomial  q,  <7(61/2)  >s  27r-pcriodic  if  and  only  if  it  can  be  written  as  an 
algebraic  polynomial  in  ei  =  e^2 ,  i.e.,  if  and  only  if  q  involves  only  even  powers,  or,  what  is  the 
same,  if  and  only  if  all  the  zeros  of  q  occur  in  symmetric  pairs.  Thus  the  quotient  r/A  in  Corollary 
5.14  can  be  equivalently  characterized  by  the  lack  of  symmetric  zeros  in  p/q. 
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If  we  take  for  <f>  a  cardinal  B-spline  and  for  77  its  2-dilate,  then  the  half-shifts  of  77  are  linearly 
independent.  In  this  case,  the  spline  wavelet  rp  of  Chui  and  Wang  (given  by  Theorem  5.5)  is 
the  minimally  supported  wavelet  for  W  guaranteed  by  Corollary  5.11  because  the  function  r  of 
Corollary  5.14  is  known  to  have  no  27r-periodic  polynomial  factor.  It  thus  follows  that  tp  has 
linearly  independent  shifts. 

6.  An  example  of  non-stationary  decompositions:  exponential  B-splines 


We  have  carried  out  the  analysis  in  this  paper  without  making  the  assumption  that  77  is  the 
2-dilate  of  <f>.  The  reason  for  this  is  two-fold:  first,  the  assumption  77  =  <p( 2-)  does  not  simplify 
either  the  idea  or  the  details  of  our  approach.  Second,  and  more  importantly,  there  are  various 
interesting  examples  where  the  “finer”  function  77  is  not  obtained  from  <p  by  dilation.  This  is  the  case 
for  example  for  exponential  B-splines,  exponential  box  splines,  and  various  radial  basis  functions. 
In  this  section,  we  briefly  discuss  what  seems  to  be  the  simplest  example  in  this  direction:  the 
exponential  B-splines. 

The  exponential  B-spline  N\  :=  N\(-\0,  ...,n)  is  a  generalization  of  the  (polynomial)  B-spline 
1V(-|0,  ...,n).  It  can  be  defined  by  its  Fourier  transform  as  follows.  Let  A  be  a  parameter  vector 
(Aj,...,An)  €  Cn.  Then 

"  eAm~'v  -  1 

ii  -  •»)  * 

The  polynomial  B-spline  corresponds  to  the  choice  A  =  0.  Splines  in  tension  correspond  to  the 
choice  n  =  4,  Ai  =  A2  =  0,  A3  =  —  A4. 

In  general,  Nx  is  (n  —  2)-times  continuously  differentiable  and  is  supported  on  [O..71].  On 
each  interval  [j . .  j  +  1],  Nx  coincides  with  a  function  in  the  kernel  K\  of  the  differential  operator 
V  :=  n"=i(£  —  Am).  The  shifts  of  Nx  are  linearly  independent  if  and  only  if 

(6.1)  Am  -  A,  i  27riS\0,  Vm,j. 


Nx(  y) 


Furthermore,  when  (6.1)  holds,  every  /  €  K\  can  be  expressed  as  a  linear  combination  of  the  shifts 
of  Nx. 

With  the  above  knowledge  in  hand,  it  should  be  clear  that  Nx  cannot  be  written  in  terms  of 
its  2-dilate,  unless  A  =  0:  upon  dilating  N\  we  obtain  a  function  which  is  piecewise  in  A^a  and 
therefore  every  element  of  «S1(ArA(2-))  is  piecewise  in  A"2a  while  Nx  is  piecewise  in  Kx.  Thus,  the 
usual  framework  of  multiresolution  analysis  cannot  be  applied  to  exponential  B-splines. 

On  the  other  hand,  from  the  Fourier  transform  of  Nx,  we  see  that 


(6.2) 


£»(») = n 

m  =  l 


g^m/2  — iy/2  j 
2 


n 


eAm/2-i'jf/2  _  j 

(Am/2  -  ty/2) 


The  second  factor  on  the  right  hand  side  of  (6.2)  is  recognized  as  IV a/2(-/2)>  and  thus 


Nx  =  Ax/2Nx/2(-/2), 
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with  A\  the  ^-periodic  trigonometric  polynomial 


A\(y)  :=  n 

m=l 


g>m-iy/2  _j_  J 
2 


Note  that  JVA/2(-/2)  is  the  Fourier  transform  of  2iVA/2(2-)  which  is  supported  on  [0..n/2],  and  is 
piecewise  in  K\  (with  breakpoints  at  the  half-integers). 

We  fix  a  vector  A  and  define  the  spaces  Sk  :=  Sk(<f>k),  k  €  TL,  with  <f>k  :=  ^a/2*( 2fc-).  The 
generators  <j>k  then  satisfy  the  non-stationary  refinement  equations 

4>k  =  2AA/2‘  +  »(72fc)0A:+l- 

We  observe  that  2Aa/2*+»(72*)  can  be  written  as 


3=0 


Ai/2* 


An/2*)e-.iy/2*+1 


where  aj(ti,...,tn)  is  the  homogeneous  symmetric  polynomial  of  degree  j  in 

The  scale  of  spaces  Sk,  k  £  7L,  clearly  satisfies  condition  (1.3)(i)  of  multiresolution.  Since 
supp^it  =  IR,  (1.3)  (ii)  follows  from  Theorem  4.3.  According  to  Theorem  4.9,  the  space  Y  :=  C\kSk 
has  dimension  <  1.  The  following  theorem,  which  is  a  special  case  of  Theorem  8.4,  provides  a 
complete  description  of  this  space. 

Theorem  0.3.  Let  {5fc}  be  a  multiscale  of  spaces  generated  by  exponential  B-splines.  Let  Y  := 
flj tez«Sfc-  Then  Y  is  one-dimensional  if  and  only  if  Re \j  0,  j  =  l,...,n.  Otherwise,  Y  is  trivial. 
In  case  Y  is  one-dimensional,  it  is  spanned  by  the  Green’s  function  G  (or  more  precisely  the 
fundamental  solution  of  the  differential  operator  V)  whose  Fourier  transform  is  given  by 


(6.4)  G(y)  =  J]  (Am  -  iy)  l. 

m= 1 


In  this  case,  Sk  =  Sk(G)  for  every  k. 

For  convenience,  we  define  from  now  on  W_0 0  :=  C\kSk,  and  obtain  in  this  fashion  the  decom¬ 
position 

I2(fil)  =  ©  Wk, 

-oo<fc<  +  oo 

valid  for  the  wavelet  decomposition  based  on  any  exponential  B-spline. 

An  interesting  and  important  problem  in  the  context  of  non-stationary  subdivisions  is  the 
stability  question.  Let  ipk  be  the  compactly  supported  wavelet  function  given  by  (the  appropriate 
scaled  version  of)  Theorem  5.5  for  4>k  :=  Arx/2‘(2fc-)-  Then  the  wavelet  space  Wk  :=  Sk+l  QSk  is  a 
principal  2-fc-shift-invariant  space  generated  by  ipk ■  The  2-fc-shifts  of  <f>k  are  linearly  independent 
ifF 

(6.5)  Am  -Xj  £2fc+17ri7Z\0, 
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as  can  be  easily  concluded  from  (6.1)  by  rescaling.  We  see  that  for  large  enough  k  (say,  k  >  ko  > 
-oo),  4>k  is  always  a  linearly  independent  generator,  and  in  particular  a  stable  generator  of  its 
space.  Further,  if  A  €  ntn,  then  the  linear  independence  (and  hence  the  stability)  holds  for  all  k. 
By  Theorem  5.5,  ipk  provides  a  stable  basis  for  the  wavelet  space  Wk  for  every  k  >  ko.  At  the  same 
time,  the  sum  ® Wk  is  orthogonal,  a  fortiori  the  sum  ©*>fco  Wk  is  orthogonal.  Nevertheless,  these 
arguments  do  not  imply  that  {ipk(-  -2-kj)}k>k0,je%d  forms  a  stable  basis  for  ®fc>fco  Wk ,  since  one 
still  needs  to  show  that  the  stability  constants  associated  with  the  basis  {ipk(-  -  2 -*j)}j€3<i  of  Wk 
can  be  chosen  independently  of  k  >  ko.  This  question  does  not  arise  in  the  stationary  case,  since 
then  tpk  is  obtained  by  dilating  V’o  and  the  stability  constants  do  not  change  with  k. 

The  main  tool  in  this  discussion  of  stability  is  the  following  consequence  of  Theorem  5.5  and 
Remark  5.8: 

Corollary  6.6.  Let  ( Sk  =  Sk(<j>k))k  be  a  nested  sequence  of  spaces  in  L 2(IR),  and,  for  k  £  TL, 
define  xpk  '=  2k/2ip(2k-),  with  rp  the  wavelet  generator  of  Theorem  5.5  corresponding  to  the  choice 
<f>  :=  <j>k(-/2k),  t)  :=  <f>k+ 1 (-/2fc)  in  that  theorem.  Let  ijk  •'=  <Pk(‘/2k),  k  £  TL,  and  — oo  <  ko  <  ki  < 
oo.  Then,  the  set 

^  '•=  {^fc(-  -  ^~kj)}k0<k<k1,je'Zd 

is  an  Li-stable  basis  for  the  space  © k0<k<kl  (with  Wk  :=  <Sfc+1  QSk)  if  and  only  if  there  exist 

positive  constants  C\  and  C2  such  that 


H^llLooCir)  —  Cj  an<i  II 1/ ^kllLooflT)  5:  1/Cl, 

for  every  k0  <  k  <  k\.  Furthermore,  the  stability  constants  Cj(9)  for  the  choice  F  :=  $  in  (1.5) 
then  can  be  taken  as  Cj,  j  =  1,2. 


Proof.  As  explained  in  the  paragraph  preceding  this  corollary,  we  only  need  to  check,  for 
each  k,  the  stability  constants  associated  with  the  basis  Vk  •—  ($*(•  —  2~fcj))j€2;j  for  Wk .  By 
Theorem  5.5  and  Remark  5.8,  these  constants  are  determined  by  the  constants  associated  with  the 
sequence  ( 4>m( •  -  2 -mi))j6Zd,  m  =  k, k  +  1.  By  scaling,  these  latter  constants  are  observed  to  be 
identical  with  the  constants  associated  with  the  sequences  (qm(-  —  m  —  k,k+l.  With  this, 

the  bounds  Cj,  j  =  1,2,  follow  from  Remark  5.8.  4 

Corollary  6.7.  Let  ( Sk(<f>k))k  be  a  nested  sequence  of  exponential  B-spline  spaces,  i.e.,  <pk  = 
^rA/2*(2fc<)  for  some  (fixed)  A  £  Cn.  Let  ko  >  -oo  be  chosen  such  that  (6.5)  holds  for  every  k  >  k o- 
Let  ip0  be  the  generator  of  Wo  defined  by  Theorem  5.5,  and  let  xpk  be  the  analogous  generator  of 
Wk,  k  £71.  Let 

®  :=  OM*  -  2  i))k>ko,j€Kd' 

Then  V  forms  a  stable  basis  for  ©fc>fco  Wk- 

Proof.  We  observe  that  rjk  '■=  <J>k(2~k-)  is  the  function  AA/2t,  and,  by  the  assumption 
here,  the  shifts  of  each  qk  form  a  stable  basis  for  5(i?fc).  All  the  functions  rjk,  k  £  TL,  are  supported 
in  [0 .  .n]  and  they  converge  uniformly  as  k  -+  oo  to  the  polynomial  B-spline  No-  From  this  it  easily 
follows  the  ?jk  converges  uniformly,  as  k  -*  oo,  to  Nq.  Thus,  for  sufficiently  large  ki,  and  for  every 
k>ku 


ll'?fc||L00(Tr‘')  <  ll^oll^en"1)  +£> 
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and 


ll1/^A:!|JL0O(Tr‘*)  ^  ll1/^o||JL00{'lT‘*)  + 

It  follows,  thus,  that  supfc>fco  ll^Jtl|£,00('ir<* ) »  and  supJk>fco  ||l/%i|t0O(Tr’1)  are  finite,  and  our  claim 
follows  from  Corollary  6.6.  6 

A  more  subtle  analysis  is  required  in  the  consideration  of  the  stability  of  the  full  basis 
—  2~kj))k  ..  We  omit  these  details  here. 

7.  Multivariate  prewavelets 

We  have  given  in  §3  various  sets  of  generators  for  the  wavelet  space  W.  In  particular,  we 
have  shown  how  to  obtain  generating  sets  which  provide  an  Testable  basis  or  more  generally  an 
orthonormal  basis  for  W.  However,  our  constructions  were  lacking  in  the  following  sense:  If  77  has 
compact  support,  then  the  elements  in  the  generating  sets  which  provide  an  Testable  basis  need 
not  be  of  compact  support,  nor  can  they  be  shown  to  decay  at  any  rate.  On  the  other  hand,  it  has 
been  proved  by  Meyer  [Me;  Chapter  III,  §6]  (and  also  Jia  and  Micchelli  [JM1])  that,  under  some 
general  assumptions  on  the  generator  <j>  of  <S°  (e.g.,  <}>  is  compactly  supported  and  provides  a  stable 
basis  for  S° ),  there  always  are  generating  sets  consisting  of  nicely  decaying  functions  which  provide 
an  I2-stable  basis  for  VV.  However,  the  proofs  of  these  facts  are  not  constructive,  hence  leave  open 
the  question  of  how  to  obtain  such  generating  sets  explicitly.  We  shall  not  provide  a  solution  to 
this  problem  in  its  entirety,  but  we  shall  build  on  previous  constructions,  of  Lorentz  and  Madych 
[LM]  and  Riemenschneider  and  Shen  [RS],  which  can  be  applied  in  certain  special  but  important 
cases. 

We  assume  throughout  this  section  that  <f>  and  77  are  I2(lR<<)-functions  that  satisfy 


(7.1) 

C  S'(r 7) 

and 

(7.2) 

supp  <t>  =  supp  77  =  IR  . 

As  before,  we  denote  the  refinement  mask  by  A,  i.e., 

4>  =  Arj. 

The  refinement  mask  plays  a  major  role  in  the  context  of  orthogonal  wavelets  (cf.  (5.9)).  But,  as 
already  observed  in  Theorem  5.5,  the  construction  of  prewavelets  is  based  on  the  function 

(7.3)  B  :=  [rj,<i>]  =  At)2, 

and  for  that  reason  we  assigned  it  the  above  special  notation,  B. 

The  derivations  of  generators  and  bases  for  W  that  were  carried  out  in  §3  involved  only  the 
function  <f> .  In  order  to  construct  stable  bases  for  W  from  (f>  that  imitated  the  decay  properties  of 
4>,  we  would  have  to  assume  that  4>  has  Testable  half-shifts,  and  this  is  a  restrictive  assumption, 
and  applies  only  to  non-stationary  refinements  (see  the  next  section).  Thus,  we  change  our  focus 
from  <p  to  r 7,  under  the  assumption  that  the  half-shifts  of  the  new  generator,  77,  are  stable.  Indeed, 
it  is  the  Z.2-stability  of  17  which  will  allow  the  construction  of  an  Testable  basis  for  W. 

We  recall  the  operator  Qo  of  (2.30). 
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Corollary  7.4.  Assume  that  <f>  and  rj  satisfy  (7.1)  and  (7.2).  A  necessary  and  sufficient  condition 
that  w  €  ^2(01^)  be  in  W  is  that  there  is  a  4 n-periodic  function  r  such  that 


(7.5) 

W  =  TTj 

and 

(7.6) 

Qo(tB)=  Y,  (tB)(>  +  v)  =  0. 

v£4irV 


If  77  has  compact  support,  then  a  sufficient  condition  that  w  have  compact  support  is  that  r  is  a 
trigonometric  polynomial  (of  period  Air).  Moreover,  this  last  property  characterizes  the  compactly 
supported  elements  ofW,  whenever  t)  has  linearly  independent  half-shifts. 

Proof.  The  first  equality  in  (7.6)  is  merely  the  definition  of  Qo .  As  for  the  second,  since 
W  C  Sl{r]),  any  function  in  W  has  Fourier  transform  of  the  form  (7.5).  Since,  for  any  functions 
f,9  €  £2(11^),  [ f,g ]  =  Qo([f,g]),  we  conclude  that 

\w,4>]  =  Qq{{tt),4>X)  =  Qo{rB). 

Since  w  €  W  if  and  only  if  [u;,  <f>\  =  0,  the  main  claim  of  this  corollary  follows. 

If  t  is  a  polynomial  and  77  is  compactly  supported,  then  rrj  certainly  is  the  Fourier  transform 
of  a  compactly  supported  function.  In  case  77  has  linearly  independent  half-shifts,  Theorem  1.3  in 
[BR]  implies  that  /  =  7-77,  with  r  a  trigonometric  polynomial,  whenever  /  €  Sl{rj)  is  of  compact 
support.  4 

With  Corollary  7.4  in  mind,  we  would  like  to  find  a  set  Vo  C  V  of  cardinality  2rf  -  1  and 
47r-periodic  functions  ru,  v  €  Vo,  that  satisfy  Qo(rvB)  =  0.  Then,  the  functions  rvrj,  v  €  Vo,  are 
in  W .  Under  certain  conditions  we  can  choose  the  r„,  v  €  Vo,  so  that  the  wv  :=  (tvtj)v,  v  €  Vo, 
provide  an  Za-stable  basis  for  W.  We  begin  by  generalizing  a  construction  used  by  Lorentz  and 
Madych  [LM]  (see  also  [JM]  and  [So]). 

We  can  decompose  the  function  B  =  [77,^]  into  its  27r-periodic  components,  as  in  (2.31): 

(7.7)  B  =  e-*Bv,  Bv  :=  QV(B)/ 2d. 

v€V 


If  77  and  <j)  are  of  compact  support,  then  (by  the  half-shift  analog  of  Lemma  2.8),  B  =  [77,  d>]  is  a 
47r-periodic  polynomial.  In  such  a  case,  the  functions  Bv  are  27r-periodic  polynomials. 

For  our  first  construction,  we  assume  that  B  is  bounded,  and  that,  for  some  uo  €  V ,  BVo  is 
bounded  away  from  zero  a.e.,  and  set  Vo  :=  V\  {uo}.  These  requirements  are  fulfilled,  for  example, 
for  Vo  =  0,  hence  Vo  =  V' ,  whenever  77  has  Z2-sta-ble  half-shifts,  and,  further,  the  47r-periodic 
refinement  mask  A  is  real,  nonnegative  and  continuous,  with  no  27r-periodic  zeros. 

Theorem  7.8.  Assume  that  <f>  and  77  satisfy  (7.1)  and  (7.2),  and  let  Bv  be  as  in  (7.7).  Let  vo  €  V, 
Vo  :=  Vr\{vo}.  Then  the  functions 

(7.9)  tv  :=  eVoBv  evBVo,  v  €  Vo, 
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satisfy  Qq{tvB)  =  0,  hence  the  functions  wv,  v  £  Vo,  with  Fourier  transform  wv  :=  rvrj,  are  in  W , 
provided  that  wv  £  If  V  and  <f>  have  compact  support,  then  the  wv,  v  £  Vo,  have  compact 

support  as  well.  If  77  has  Li-stable  half-shifts,  and  both  B  and  1  /BVo  are  essentially  bounded  on 
TT11,  then  (wv)veVo  provides  an  Li-stable  basis  for  W. 

Proof.  Let  v,u£V.  Since  Bv  is  27r-periodic, 


Qo(euBvB)  =  BvQo(euB)  =  BVQU(B )  =  BvBu2d. 

Application  of  this  equality,  once  with  v  =  v,  u  =  vq  and  then  with  the  opposite  choice,  proves 
that 

2~dQ0(TvB )  =  BVBV0  -  BV0BV  =  0. 

Hence  (7.6)  is  satisfied  and  the  functions  wv,  v  £  Vo,  are  in  W. 

If  rj  and  <j>  have  compact  support,  then,  by  the  half-shift  analog  of  Lemma  2.8,  B  =  {77,  </>]  is  a 
trigonometric  polynomial,  hence  so  is  each  Bv  and  each  r„.  This  implies  that  each  wv  is  compactly 
supported. 

To  show  the  /^-stability  of  the  (iVv)v€V,  we  consider  the  matrix  T  :=  (ru,u)vu6U  with 
diagonal  elements  rUt„  :=  —B „0,  v  £  V,  and  with  off-diagonal  elements  rV  Vo  =  Bv,  v  £  Vo,  and 
with  all  other  entries  zero.  We  observe  that  T [cvfj)v&v  coincides  with  {wv)veV°  in  all  the  Vo- 
entries,  and  therefore,  for  proving  the  desired  stability  it  suffices  to  show  that  the  shifts  of  the 
inverse  transforms  of  T(e„^)vGV  are  stable.  Recall  that  we  are  assuming  that  the  half-shifts  of  17 
are  stable,  or,  equivalently,  that  the  full-shifts  of  (??(•  +  V))vey  are  stable.  Thus,  by  (iv)  of  Theorem 
2.26,  it  remains  to  show  that  [|T||  and  ||T_,(|  are  essentially  bounded  on  TTd.  Since  we  assume 
that  B  is  bounded,  so  is  each  component  Bv,  hence  T  has  all  entries  bounded.  On  the  other  hand, 
|  det  T|  =  |RUo|2*  hence  by  our  assumption,  is  bounded  away  from  0.  This  implies  that  both  ||T|| 
and  ||T-1||  are  bounded  a.e.  4 

We  note  that  the  boundedness  assumption  on  B  is  automatically  satisfied  whenever  the  full- 

shifts  of  77  and  <p  are  stable,  since  \B\  =  rj<t>,  with  each  factor  on  the  right  being  bounded  because 
of  the  stability  assumption. 

The  other  assumptions  of  Theorem  7.8  are  also  met  in  many  instances.  The  most  important 
example  is  recorded  in  the  following  corollary,  which  also  admits  straightforward  extensions  to  the 
non-compact  support  case. 

Corollary  7.10.  Let  <f>  and  77  be  two  compactly  supported  functions  with  the  half-shifts  of  77 
and  the  full-shifts  of  4>  being  Li-stable.  Suppose  that  (7.1)  and  (7.2)  are  satisfied,  and  that  the 
refinement  mask  A  is  a  trigonometric  polynomial  (continuity  of  A  would  suffice,  as  well).  Then 
4>  1  :=  4>  *  </*(-•)  rt.nd  771  :=  77  *  77(--)  satisfy  the  conditions  of  Theorem  7.8  for  Vo  =  0,  and  hence  the 
sequence  (wv)vey,  of  (compactly  supported )  functions  defined  there  (with  respect  to  <f>\  and  771) 
forms  a  stable  basis  for  the  corresponding  wavelet  space. 

Proof.  We  observe  that  771  =  |^j2  and  <P  1  =  |<£|2>  and  hence  <f>\  —  Because  of  the 

compact  support  of  r]\  and  <f>i,  the  non-negative  function 

B  :=  [T?i,d>iJ  =  |A|2  77j2 
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is  a  trigonometric  polynomial  (Lemma  2.8),  hence  bounded  (alternatively,  it  is  bounded  because 
of  the  stability  assumptions  on  rj  and  <f>).  Since  the  half-shifts  of  rj  are  stable,  rj  possesses  no  4x- 
periodic  zero,  hence  neither  does  rjx  =  |tj|2.  Consequently,  rjl  vanishes  nowhere.  Also,  because  of 
the  stability  of  the  shifts  of  <f>,  A  has  no  2^-periodic  zeros  (since  such  zeros  would  be  inherited  by 
(f>,  hence  by  <f>).  This  means  that  B  is  a  non-negative  47r-periodic  function  without  any  27r-periodic 
zeros.  Consequently,  Bq  is  a  strictly  positive  (2;r-periodic)  trigonometric  polynomial.  Now  apply 
Theorem  7.8.  ♦ 

We  next  describe  a  general  procedure  for  constructing  functions  r  which  satisfy  (7.6).  The 
vertices  V  form  a  group  under  addition  modulo  one.  If  J  is  one  of  its  subgroups,  then  the  distinct 
cosets  v  +  J  form  a  partition  of  V .  We  let  R  C  V  be  a  set  of  representers  for  these  distinct  cosets. 
A  partition  R  —  l?oUJ?i  into  disjoint  sets  gives  the  sets  Kj  :=  Uugj?y(v  +  J),  j  =  0, 1,  which  are 
a  partition  of  V.  Note  that  if  ev,  v  €  V,  is  an  exponential  which  is  not  constant  on  4irJ,  then 

£„€4*  =  0. 

Theorem  7.11.  Assume  that  <}>  and  tj  satisfy  (7.1)  and  (7.2),  and  define,  as  before,  B  := 

Let  J  be  any  subgroup  of  the  group  V,  let  v  be  any  element  ofV  for  which  ev  is  non-constant  on 
4 ttJ,  and  let  K  be  any  union  of  cosets  (in  V)  of  J  which  contains  0.  Then  the  function  wv  k,  w>th 
Fourier  transform 

wv,i<  :=  evrj  1 ?(•  +  a), 

a€4irA'\0 

is  in  W  provided  it  is  in  L2(iRd).  Moreover,  if  ij  and  <f>  have  compact  support,  then  w  is  also  of 
compact  support. 

Proof.  In  case  tj  and  <f>  are  compactly  supported,  B  =  [77,  is  a  trigonometric  polynomial, 
and  hence  each  inVi/<  (which  is  then  a  well-defined  /^-function)  is  compactly  supported. 

To  prove  the  main  claim  of  this  theorem,  it  is  enough,  in  view  of  Corollary  7.4,  to  show  that 

r:=  ]~[  B(-  +  a) 

a£4irK 


satisfies 

^2  ev{-  +  p)r{-  +  p)  =  0. 

Since  r(-  +  1/)  =  r,  v  6  47r  J,  we  can  write  this  last  sum  as 

Y  ett(-  +  r  +  i/)r(-  +  r  +  i/)  =  e«(’  +  rM*  +  r)  €v^' 

r £4nfi  t/£4irj  r£4r/i  i/£4rj 


The  last  sum  is  0  because  ev  is  not  constant  on  4nJ.  4 

The  choice  K  —  V  in  the  last  theorem  shows  that  all  the  (Zd/2)\Zd-translates  of  the  function 
in,  defined  by 

w  :=  rj  fj  B(-  +  a), 

a€4»rV' 
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are  in  W ,  provided  that  w  €E  Li-  It  is  then  easy  to  prove  that  the  F'-shifts  of  w  provide  a  basis 
for  the  wavelet  space.  There  is  a  close  relation  between  the  function  w  here  and  the  generator 

/N 

w  of  W  of  Corollary  3.19,  only  that  there  we  used  <f>2  =  [4>,<t>],  while  here  we  use  the  function 
B  =  It  follows,  for  example,  that  if  <f>  and  77  are  compactly  supported  and  the  refinement 

mask  A  is  a  polynomial,  w  here  enjoys  a  smaller  support  than  w  of  Corollary  3.19.  However,  unless 
the  half-shifts  of  <f>  are  stable,  neither  of  these  generators  is  expected  to  provide  a  stable  basis  for 
W. 

The  simplest  instance  of  Theorem  7.11  occurs  when  J  :=  {0,o}  is  a  group  of  order  2.  We 
obtain  the  following  extension  of  Theorem  5.5  to  the  multivariate  setting.  Here,  as  before,  B  is 
defined  as  [rj,  4>\  —  A j?2- 

Corollary  7.12.  Assume  that  (f>  and  77  satisfy  (7.1)  and  (7.2).  If  v  £  V'  and  a  €  47rV"  satisfy 
e„(a)  =  -1,  then  the  function  w  with  Fourier  transform 


w  =  ev  B(-  +  a)rj 


is  in  W  provided  it  is  in  /^(IR^).  Moreover,  if  77  and  <t>  are  of  compact  support,  then  w  is  also  of 
compact  support. 

In  some  instances,  it  is  possible  to  find  an  Z,2-stable  basis  from  among  the  functions  of  Corollary 
7.12,  as  is  shown  in  the  following  theorem  of  Riemenschneider  and  Shen  [RS]  (see  also  [JM],  [RSI] 
and  [CSW]): 

Theorem  7.13.  Assume  that  <f>  and  rj  satisfy  (7.1)  and  (7.2)  and  that  B  :=  (fj,  4>\  is  real-valued. 
Assume  that  4>  has  L2 -stable  full-shifts  and  ri  has  Lo-stable  half-shifts.  Assume  further  that  there 
is  a  one-one  mapping  a  from  V'  to  47rV'  that  satisfies  the  following  two  conditions: 

(a)  eu(a(u))  =  -1,  for  every  v  £  V'; 

(b)  ev_u(a(u)  -  a(u))  =  -1,  for  all  v,  u  €  V',  unless  v  =  u. 

Then,  the  functions  wv,  v  €  V' ,  defined  by  their  Fourier  transforms 

wv  :=  2dev  B( •  +  Q(v))rj,  v  €  V\ 

provide  an  L2 -stable  basis  for  W.  Furthermore,  if  the  full-shifts  of  <p  are  orthonormal  and  the 
half-shifts  of  77  are  also  orthonormal,  then  (wv)v£V,  provides  an  orthonormal  basis  for  W.  If  rj  and 
<j>  have  compact  support,  then  the  functions  wv,  v  €  V' ,  are  also  of  compact  support. 

Proof.  It  is  easy  to  conclude  from  the  stability  assumption  on  77  that  each  wv,  v  £  V' ,  is 
in  L 2(1R'1),  and  hence  each  wv  is  a  well-defined  I2-function.  From  Corollary  7.12  and  assumption 
(a),  we  conclude  that  each  wv  is  in  W.  This  corollary  also  implies  that  wv  is  compactly  supported 
whenever  77  and  <f>  are. 

We  introduce  the  functions  te*,  v  £  V'  with  Fourier  transform  iD*  =  wv/ri.  These  functions 
are  in  because  rj  is  bounded  away  from  zero,  thanks  to  the  stability  assumption  on  the 

half-shifts  of  77.  We  now  compute  the  Gramian  of  these  functions.  First,  we  see  that 

2~u[w'v,  tt£J  =  +  a(v))B(-  +  a(7r))[7,  =  ev_uB(-  +  a(r))r(-  +  a(u)). 
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(Here,  we  wrote  B( ■  +  a(u))  is  instead  of  B(-  +  a(u)),  since  B  is  assumed  to  be  real.)  Therefore, 


(7.14)  [wl,K}=  Y  ^uK*  + t2)  =  22d  Y  ev-u(- +  +  <*(v)  +  v)B(- +  a(u)  +  fi). 

#j€4irV  #i€4jrVr 


For  any  /*  £  4ttV,  the  terms  in  (7.14)  corresponding  to  (i  and  pi  +  a(v)  -  a(u )  are  negatives 
of  one  another  because  of  our  assumption  (b),  and  the  (a(v)  -  a(u))-periodicity  of  the  term  B( •  + 
a(u)  +  n)B{-  +  a(u)  +  /r).  Hence  [u;*,  t£*]  =  0,  v  ^  u.  On  the  other  hand,  for  v  —  u, 

(7.15)  [wl,wl\  =  22d  Yj  B(.  +  a(v)  +  ti)2  =  22d  Y  B(- +  ^)2  =  22dQ0(B2). 

ti£4irV  /i£4irVr 

Since  B 2  =  A 2rji  =  |>1|2»74  =  4>2rj2  and  Qo(4>2)  =  4>2,  we  have 

Qo(B2)(x)  £  [m  . .  M]{4>{x)2  a.e., 

with  m  and  M  the  essential  infimum,  respectively,  supremum  of  rj2.  Since  both  m  and  M  are 
positive  and  finite  by  the  stability  assumption  on  tj,  while  <t>2  is  essentially  bounded  away  from  0 
and  infinity  by  the  stability  assumption  on  <j>,  we  conclude  that  also  Q0(B2)  is  essentially  bounded 
away  from  0  and  infinity.  We  thus  conclude  that  the  Gramian  associated  with  (wC)u6V"  *s  d*ag°naK 
with  the  diagonal  entries  bounded  above  and  below  by  positive  constants.  On  the  other  hand, 
VV  :=  (wv)veV  is  obtained  from  (u>l)veV,  by  multiplying  by  the  scalar  matrix  T  :=  rj2I.  Again, 

the  stability  assumption  on  the  half-shifts  of  77  implies  that  rj  and  1  /rj  are  bounded,  hence  that  ||T|| 
as  well  as  ||T—1 1|  are  bounded.  Thus,  from  Theorem  2.26(iv),  we  conclude  that  the  basis  (wv)v^v, 
is  stable. 

Finally,  when  77  has  orthonormal  half-shifts,  77  =  2~d I2  a.e.,  and  hence 


?dQo(B2)  =  2UQ0(1>2^2)  =  2dQ0(52)  =  2^2. 


If  also  4>  has  orthonormal  full-shifts,  4>  =  1  and  hence  2 2dQo(B2)  =  2d.  Thus,  (7.15)  implies  that 
[«;*,«;*]  =  2rf,  hence  [«/„,«>„]  =  ^2[u)*,5?*j  =  2— d2d  =  1,  and  we  conclude  that  G(\V)  =  /,  or 
equivalently,  that  W  is  an  orthonormal  basis.  4 

We  make  the  following  additional  remarks  concerning  Theorem  7.13.  As  Riemenschneider  and 
Shen  [RS]  note,  it  is  easy  to  construct  mappings  with  properties  (a)  and  (b)  in  the  case  d  =  1,2,3. 
However,  Riemenschneider  and  Shen  also  note  that  there  are  no  such  mappings  when  d  >  3.  On  the 
other  hand,  there  is  some  hope  that  turning  to  the  more  general  elements  of  W  given  in  Theorem 
7.11,  an  analogue  of  Theorem  7.13  may  be  established  in  higher  dimensions. 

We  have  assumed  in  Theorem  7.13  that  the  function  B  =  [77,^]  is  real.  Since  also  B  =  At)2,  B 
is  real  if  and  only  if  the  mask  A  is  real.  This  is  true  for  example  if  4>  is  real- valued  and  symmetric 
about  the  origin  and  77  =  4>(2  ).  Moreover,  the  assumption  that  B  (or  A)  is  real  can  be  somewhat 
weakened.  For  example,  the  proof  given  above  supports  the  following  claim. 
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Remark  7.16.  The  construction  detailed  in  Theorem  7.13  remains  valid  in  case  B  =  ejB'  for 
some  real  B'  and  some  j  G  7Zd/2. 


8.  Box  Splines 

Box  splines  were  introduced  by  the  first  two  authors  in  [BD]  and  their  exponential  general¬ 
ization  (sometimes  referred  to  as  “exponential  box  splines”)  was  introduced  by  the  third  author 
in  [Rl].  Box  splines  have  become  a  main  theme  in  Multivariate  Spline  Theory,  and  it  is  certainly 
beyond  the  scope  of  this  section  to  provide  a  good  account  on  box  splines.  We  do  not  even  attempt 
to  provide  an  overview  of  box  splines  in  the  context  of  wavelet  decompositions,  because  of  the 
already  rich  literature  on  that  matter.  Thus,  our  only  aim  here  is  to  illustrate  the  material  detailed 
in  previous  sections  via  a  discussion  of  this  class  of  examples. 

To  define  a  box  spline,  we  let  T  be  a  finite  index  set  consisting  of  pairs  of  the  form 

7  =  (z-pA-,),  x7  G  ffi.d\0,  €  C. 

The  box  spline  M  :=  M\  can  then  be  defined  via  its  Fourier  transform  as 

(si)  = n 

The  notation  is  indicative  of  the  fact  the  we  usually  hold  the  directions  (x7)7er  fixed,  but  may 
vary  the  parameters  A  :=  (A7)^ep.  Assuming  that  span  (z^)^gp  =  1R  (as  we  do  throughout),  the 
box  spline  is  a  compactly  supported  piecewise-exponential-polynomial  function  supported  in  the 
zonotope 

Zr  :=  {^T^f-yZ-),  :  t~,  G  [0..1]}- 

7er 

The  polynomial  case  corresponds  to  the  choice  A  =  0.  Exponential  B-splines  are  obtained  when 
d  =  1  and  z-,  =  1,  all  7.  Tensor  splines  are  obtained  whenever  all  the  directions  are  standard  unit 
vectors.  The  box  spline  is  positive  in  the  interior  of  Zp  whenever  A  is  real-valued. 

We  first  observe  that 

M\(y)  =  ^A/i(y/2)A/A/2(y/2), 


with 


Axiy)  ■•=  n 

•rer 


This  suggests  the  choice  ^ 

£*:=2-wMA/2*(-/2fc), 

since  then  4>k  -  2dAA/2*+.  (•/2*+,)0*+i.  To  ensure  the  fact  that  A*  is  a  2?r-periodic  polynomial, 
we  assume  that 


(8.2) 


z7  €  7Zd\0,  V7€l\ 
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Assuming  (8.2),  we  can  define  the  multiscale  generated  by  the  box  spline  M\  as 
(8.3)  Sk  :=  Sk(<t>k ),  4>k  “  MA/2*  (2fc-). 

As  before,  we  use  special  notation  for  fo  and  <f>\ : 

<t>:=4>0,  T}  :=  <f>\. 

Note  that  this  is  a  stationary  multiscale  if  and  only  if  A  =  0,  i.e.,  if  and  only  if  M\  is  a  polynomial 
box  spline. 

Since  each  4>k  is  compactly  supported,  we  know  that  property  (1.3)(ii)  is  satisfied  here.  With 
regard  to  (1.3)(iii),  we  have  the  following  extension  of  Theorem  6.3: 

Theorem  8.4.  Let  {<Sfc}  be  a  multiscale  of  spaces  generated  by  the  box  spline  M\  (as  in  (8.3)). 
Let  Y  :=  n*g2Z<Sfc.  Then  Y  is  one-dimensional  if  and  only  if  Re  Xy  ^  0,  *  •_  r.  Otherwise,  Y  is 
trivial.  In  case  Y  is  one-dimensional,  it  is  spanned  by  the  fundamental  solution  of  the  differential 
operator  V  :=  —  ^-y)  (where  Dx  is  the  directional  derivative  in  the  x-direction)  whose 

Fourier  transform  is  given  by 

(8.5)  G(y )  =  fJ(A^  -  ty  •  *,)“l. 

•yer 

In  this  case,  Sk  =  Sk(G),  for  every  k. 

Proof.  Let  €  C\kSk  be  a  non-zero  function.  Since  ^  €  S~k,  it  is  a  linear  combination 

of  the  2fc-shifts  of  the  box  spline  <f>k  :=  M2*A(-/ 2k).  Since  the  ratio  4>k[G  is  (a  trigonometric 
polynomial)  of  period  2fc+17r,  (2.1)  implies  that  every  function  in  Sk  must  have  the  form 

f  =  rkG, 

with  Tk  2fc+17r-periodic.  From  the  fact  that  supp  G  =  IR*,  we  conclude  that  all  rk  agree  a.e.  with  one 
measurable  function  r,  and  this  function  is  necessarily  invariant  under  all  27r-dyadic  shifts.  Lemma 
4.6  then  implies  that  r  =  const,  hence  the  Fourier  transform  of  every  function  in  the  intersection 
is  a  scalar  multiple  of  G.  Therefore,  this  intersection  is  trivial  if  and  only  if  it  does  not  contain 
G,  and  otherwise  it  is  spanned  by  G.  Since  the  ratio  <f>k/G  is  2fc+17r,  then,  again  by  (2.1),  G  €  Sk 
if  and  only  if  G  €  /^(HT*)-  Consequently,  the  proof  of  the  theorem  is  reduced  to  the  proof  of  the 
following  claim: 

UG  €  L2(IRrf)  if  and  only  if  ReA7  ^  0  for  every  7  €  T.” 

If  Re \j  0  for  every  characteristic  value  Aj,  then  we  easily  verify  that,  because  Xr  is  of  rank 
d,  G  is  in  T2,  hence  so  is  G.  On  the  other  hand,  if,  for  some  7,  ReA7  =  0,  then  G  cannot  lie  in 
X2(IR<<),  since  it  is  not  even  in  £-2(0)  whenever  the  open  set  f l  contains  points  from  the  zero  set  of 
yv-A7-t'y.x7.  ♦ 

Assuming  (8.2),  the  shifts  of  M\  are  linearly  independent  only  if  T  is  unimodular,  which 
means,  by  definition,  that  every  dx  d  matrix  whose  rows  are  taken  from  the  multiset  (x7)^€f  has 
determinant  -1,0  or  1.  Further,  if  A  is  real-valued,  the  unimodularity  assumption  is  also  sufficient 
for  linear  independence.  For  these  reasons,  we  assume  for  the  remainder  of  this  section  that  T  is 
unimodular  and  A  is  real-valued. 


37 


We  want  now  to  consider  the  possible  applications  of  the  constructions  proposed  in  the  last 
section  to  box  splines.  It  is  hard  to  apply  Theorem  7.8  directly,  since  it  requires  information  on  the 
function  B ,  while  the  available  information  here  is  on  the  mask  A.  Nevertheless,  if  A  has  the  form 


A  =  const\e_jA', 

where  j  €  TLd /2  and  A'  is  non-negative,  then  we  might  choose  v  €  V  such  that  j  -  v  €  7Ld  to 
obtain 

Bv  :=  Qv(B)/2d  =  const  ^  e„(-  +  p)e.j(-  +  p)B'(-  +  p.)  =  ev.jQ0(B'), 

where  B'  :=  (1/const )ejB.  Since  B'  =  A'rj 2,  and  A'  is  non-negative,  so  is  B' .  This,  together  with 
the  stability  assumption  on  the  half-shifts  of  q  and  the  shifts  of  4>,  implies  that  Q0(B')  does  not 
vanish,  hence  l/Bv  is  bounded,  and  we  arrive  at  the  following  conclusion: 

Corollary  8.8.  Let  M  be  a  box  spline  defined  by  a  unimodular  T  with  real  parameter  vector  A. 
Assume  that  T  satisfies  also  the  following  “parity”  condition:  “T  can  be  partitioned  into  pairs  such 
that  each  pair  (7,70  satisfies 

(x-pA-O  =  £(7>7  )(xy,— Ay), 

where  £(7,70  €  {±1}.”  Let  B  :=  [^i,d>o]>  with  <f>k  defined  as  in  (8.3).  Then  BVo  vanishes  nowhere 
on  TTd,  where  v0  £  V  is  determined  by  the  condition 

v0  =  Yl  xi/2,  mod  TLd . 

-rer 


Consequently,  the  construction  detailed  in  Theorem  7.8  can  be  applied  with  respect  to  this  r0- 

Proof.  Since  B  here  is  a  polynomial,  it  is  clear  that  the  functions  wv  defined  in  Theorem 
7.8  are  in  L-i.  Also,  for  the  same  reason,  B  is  bounded.  Thus,  to  apply  Theorem  7.8,  we,  indeed, 
need  only  to  prove  the  boundedness  of  1  /BVo.  In  view  of  the  remarks  preceding  this  corollary,  it 
suffices  to  show  that  the  mask  A  in  the  equation  d>o  =  2dA<f>\  is  of  the  form 

A  =  ejA', 

with  A'  non-negative  and  j  —  v o  €  7Zd.  Here, 


My)  =  II 

-rer 


e>,/2-ir,  y/2  y  j 
2 


Let  (7,70  be  a  pair  in  the  partitioning  of  T.  Then 

g.X.,/2  — «xvy/2  _|_  j  eA7./2-tt,.-y/2  ^  j 


(l— e^/2  +  1  +  v/*yc0sh(\^2)  +  cos(x7  •  y/2)). 
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The  second  factor  above  is  non-negative.  Multiplying  the  first  factor  over  all  pairs  (7,7')!  we  obtain 
an  expression  of  the  form 

const  aCj/2, 

where  j  :=  53e(-y,7*)=i  x~r  =  1/2  £76 r  z7  6  TLd ,  since  each  direction  appears  either  as  an  x7  or  an 
x7<  and  x7<  =  x7  if  £(7,7')  =  1  while  xy  =  -x7  otherwise.  ^ 

It  should  be  clear  that,  under  the  assumptions  of  the  last  corollary,  Theorem  7.8  can  be  applied 
to  obtain  stable  bases  for  all  the  wavelet  spaces  of  the  multiscale  generated  by  the  box  spline  M. 
Also,  the  assumption  that  A  is  real  is  convenient  but  not  essential.  In  general,  to  obtain  a  box  spline 
that  satisfies  the  above  assumptions,  one  can  start  with  any  M  that  is  defined  by  a  unimodular 
T,  and  replace  M  by  M  *  A/(  — •).  The  box  spline  obtained  in  this  way  corresponds  to  the  choice 
uo  =  0  in  the  above  corollary.  The  other  variants  can  be  obtained  the  shifting  that  box  spline  by 
j  €  V. 

If  d  <  3,  we  can  also  try  to  employ  the  construction  detailed  in  Theorem  7.13.  Here,  given  a 
unimodular  T,  we  want  the  mask  A  to  be  of  the  form 

A  =  tjA! 

for  some  real  A'.  In  the  polynomial  case  (i.e.,  when  A  =  0),  this  assumption  is  always  satisfied 
since  then  for  A  :=  A0  we  have 


My)  =  cj(y)  J!  ccs(y  •*-»)> 

-rer 


with  j  :=  -l/2^7erI7-  This  observation  immediately  extends  to  the  case  when  A  e  ilR^,  but, 
however,  does  not  extend  to  an  arbitrary  A.  On  the  other  hand,  if  M  is  a  box  spline  as  in  Corollary 
8.6,  and  A/'  is  a  polynomial  box  spline  (with  a  unimodular  set  of  direction),  then  A/  *  A/'  satisfies 
the  requirements  of  Theorem  7.13. 


We  mentioned  previously  that  for  non-stationary  subdivisions  one  should  not  exclude  the 
possibility  that  the  half-shifts  of  <f>  are  stable.  Box  splines  provide  an  excellent  illustration  of  this 
point.  In  order  to  check  the  stability  of  the  half-shifts  of  <p,  we  consider,  as  before,  the  function 


<{>  =  |Aa/2(  /2)|^-  By  our  assumptions,  the  half-shifts  of  77  are  linearly  independent,  hence  stable, 
which  means  that  rj  is  positive  on  HC*.  Therefore,  the  zeros  of  <fi  are  identical  with  those  of  the 


mask 


Ax/i(y/2)  = 

-rer 


(eXT/2-iy/2  r1  y  j) 


We  observe  that  the  factor  eXj/2  ,s'/2  x->  -j-  1  has  zeros  in  IR^  if  and  only  if  A7  =  0  (recall  that  we 
are  already  assuming  that  A7  is  real).  Thus  we  obtain  the  following  interesting  result: 


Corollary  8.7.  Let  A/*  be  the  box  spline  given  by  a  unimodular  F  and  a  real  A.  Then  M \  vanishes 
nowhere  if  and  only  if  A  contains  no  zero  entry.  Consequently,  this  last  condition  is  equivalent  to 
the  stability  of  the  half-shifts  of  M\. 
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Stronger  results  n  be  obtained  by  a  finer  analysis.  It  can  be  shown  that,  assuming  only  (8.2) 
(which  is  imbedded  in  the  last  corollary  in  the  unimodularity  assumption  on  T),  the  stability  of 
the  half-shifts  of  M\  is  equivalent  to  the  existence  of  non-trivial  functions  in  the  intersection  C\kSk, 
with  (Sk) k  the  multiscale  generated  by  the  box  spline  M\. 

The  stability  of  the  half-shifts  of  M\  leads  to  painless  constructions  of  compactly  supported 
stable  bases  for  the  wavelet  space.  Here  is  a  sample  result  in  this  direction: 

Proposition  8.8.  Assume  that  <f>  and  tj  satisfy  (7.1)  and  (7.2),  and  assume  that  the  half-shifts  of 
4>  are  stable.  Let  w  be  either  the  generator  for  the  wavelet  space  introduced  in  Corollary  3.19  or  in 
the  paragraph  after  Theorem  7.11.  Then  w  provides  (i.e.,  the  (ZZd/2)\2Zd -translates  of  w  form)  a 
stable  basis  for  W . 

The  last  result  is  less  impressive  than  it  might  look  at  first.  Indeed,  considering  the  box 
spline  multiscale  and  assuming,  say,  that  A  is  real  and  contains  non-zero  entries,  we  can  easily 
find  single  compactly  supported  stable  generators  to  each  of  the  wavelet  spaces  associated  with  the 
multiscale  {<t>k)  k  generated  by  M.  Still,  as  mentioned  already  in  section  6,  it  is  crucially  important 
to  understand  the  behavior  of  the  stability  constants  as  k  varies,  and  in  the  case  of  box  splines 
these  constants  deteriorate  fast  as  k  increases.  This  can  be  observed  as  follows:  if  we  rescale  each 
4>k  and  4>k+ 1  by  2k+l ,  and  denote  the  functions  obtained  by  r)k  and  T)k+\  respectively,  we  obtain  a 
refinement  equation  of  the  form 

rjk  =  2dA\/2k+1  tyt+i- 

Thus  complex  zeros  of  the  fcth  order  mask  converge  (exponentially)  to  the  real  domain,  as  k 
increases.  One  might  attempt  to  choose  very  large  initial  entries  for  A,  yet  the  results  of  [DR] 
indicate  that  the  asymptotic  approximation  properties  of  Sk(4>k )  deteriorate  exponentially  with 
the  growth  of  A. 

The  above  discussion  demonstrates  the  difficulty  of  controlling  the  stability  constants  in  case 
the  wavelet  constructions  are  based  on  the  stability  of  the  half-shifts  of  M\.  On  the  other  hand,  the 
constructions  that  make  use  of  the  refinement  equation  (such  as  the  one  detailed  in  Corollary  8.6) 
require  only  the  stability  of  the  half-shifts  of  2(2-).  Using  methods  similar  to  those  employed 
in  §6,  it  can  be  shown  that  for  such  constructs  the  stability  constants  do  not  blow  up  as  k  — ►  oo. 
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